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ABSTRACT 


This  system  study  extrapolates  present  physics  knowledge  and  technology  to  predict 
the  anticipated  characteristics  of  D-®He  spherical  torus  fusion  reactors  and  their 
sensitivity  to  uncertainties  in  important  parameters.  Reference  cases  for  steady- 
state  1000  MWe  reactors  operating  in  H-mode  in  both  the  1st  stability  regime  and 
2nd  stability  regime  were  developed  and  assessed  quantitatively.  These  devices 
would  have  a  very  small  aspect  ratio  (A=1.2),  a  major  radius  of  about  2.0  m,  an  on- 
axis  magnetic  field  less  than  2  T,  a  large  plasma  current  (80-120  MA)  dominated  by 
the  bootstrap  effect,  and  high  plasma  beta  (>  0.6).  The  estimated  cost  of  electricity 
is  in  the  range  of  60-90  mills/kW-hr,  assuming  the  use  of  a  direct  energy  conversion 
system.  The  inherent  safety  and  environmental  advantages  of  D-®He  fusion  indicate 
that  this  reactor  concept  could  be  competitive  with  advanced  fission  breeder  reactors 
and  large-scale  solar  electric  plants  by  the  end  of  the  21st  century  if  research  and 
development  can  produce  the  anticipated  physics  and  technology  advances. 
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CHAPTER  1 

Introduction  and  Historical  Review 
1.1  Thesis  Objectives 

This  study  investigates  the  potential  of  a  spherical  toms  fusion  reactor  fueled 
with  Deuterium  and  Hehum-3,  isotopes  of  hydrogen  and  hehum,  for  commercial  elec¬ 
tric  power  production.  The  reactor  system  study  contained  in  this  thesis  presupposes 
that  the  present  knowledge  of  plasma  physics  and  nuclear  fusion  technology  may  be 
extrapolated  to  predict  the  anticipated  characteristics  of  future  reactors.  The  result¬ 
ing  conclusions  must  be  considered  speculative  because  numerous  gaps  in  the  physics 
knowledge  and  significant  technologv  m  uncertainties  still  exist.  Nevertheless,  the 
results  of  this  study  are  extremely  useful  because  they  identify  critical  parameters 
or  assumptions,  show  which  are  important  physical  or  technological  constraints,  and 
indicate  the  broad  directions  in  which  future  research  and  development  should  pro¬ 
ceed.  The  inclusion  of  unit  costs  allows  a  qu^ultitative  estimate  to  be  made  of  the 
economic  potential  of  this  fusion  reactor  concept  relative  to  other  energy  systems. 

The  reactor  system  study  presented  herein  uses  parametric  investigations  to 
optimize  the  reactor  configuration  by  allowing  a  wide  range  of  physical  parsuneters 
to  vary  while  bounded  by  physics  and  engineering  constraints.  Besides  showing 
sensitivities  associated  with  specific  parameters,  these  studies  axe  valuable  in  three 
ways:  (1)  they  evaluate  alternative  confinement  geometries  eis  potentied  reactor  sys¬ 
tems,  thus  aiding  the  planning  of  more  effective  plasma  physics  research  programs; 
(2)  they  evaluate  the  technological  and  engineering  problems  involved  with  the  sys¬ 
tem,  thus  aiding  the  planning  of  more  effective  fusion  technology  research  programs; 
and,  (3)  they  evaluate  the  operational,  economic,  and  environmental  advantages  of 


1 


2 


fusion  energy,  thus  providing  justification  for  more  effective  fusion  research  and  de¬ 
velopment  programs.  An  important  objective  of  this  reactor  study  is  therefore  to 
quantitatively  assess  the  potential  of  a  D-®He  spherical  torus  fusion  reactor. 

1.2  D-^He  Fusion 

Two  major  motivations  for  developing  fusion  energy  are  the  potential  to  be 
safer  and  more  environmentally  benign  than  the  competition  and  the  potential  to 
secure  a  virtually  unlimited  future  energy  source.  For  the  long  term  (i.e.,  22nd 
century  and  beyond),  fusion  energy  systems  must  demonstrate  attractive  economic, 
safety  and  environmental  chciracteristics  comparable  with  advanced  fission  breeder 
reactors  and  large-scale  solar  electric  plants  [1,  2].  Ideally,  a  purely  aneutronic 
fusion  reactor  based  on  the  exotic  ^He-®He  fuel  cycle  may  one  day  prove  to  be 
viable.  At  present,  however,  the  most  promising  candidate  for  an  energy  source 
is  a  neatly  aneutronic  D-®He  fusion  reactor.  Numerous  studies  of  advanced  fuel 
cycles  and  fusion  reactors  have  been  completed  [3,  4,  5]  and  the  current  ARIES-III 
study  is  conducting  a  detailed  analysis  of  a  conceptual  D-®He  fusion  reactor  based  on 
advanced  tokamak  parameters  [6].  Compeued  with  advanced  fission  breeder  and  D-T 
fusion  reactors,  a  D-^He  reactor  offers  substantial  advantages  which  could  eventually 
prove  to  be  competitive. 

The  most  significant  advantages  of  a  D-®He  fuel  cycle  are  reduced  neutron 
activation  and  tritium  handling  problems.  The  basic  fusion  reaction  involved  in  this 
fuel  cycle  is 

P  -f-3  ^  +  p  -h  IS.ZMeV 

The  important  points  to  note  are  that  both  species  involved  in  this  reaction  are 
nonradioactive  and  that  100%  of  the  energy  released  in  this  reaction  is  in  the  form 
of  nonradioactive  charged  particles  which  deposit  their  energy  in  the  plasma  and  help 
maiintain  the  plasma  energy  balance.  This  fuel  cycle  is  conditionally  aneutronic  for 
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ion  temperatures  around  50  keV  and  ^He;D  fuel  ratios  greater  than  4.0,  but  this  ratio 
may  be  too  high  for  plasma  ignition  and  operation.  For  a  more  realistic  fuel  ratio  of 
about  1:1,  some  neutrons  will  be  produced  directly  by  the  D(d,n)®He  reaction  and 
indirectly  by  the  D(d,p)T,T(d,n)^He  and  D(d,p)T,T(t,2n)  ^He  reaction  chains,  but 
these  account  for  less  than  6%  of  the  total  fusion  power  compared  with  about  80% 
for  the  D-T  fuel  cycle.  The  reduced  neutron  irradiation  of  the  reactor  first  wall  and 
shield  will  result  in  longer  component  service  life  and  lower  activation  levels,  thus 
minimizing  maintenance  problems  and  shielding  requirements.  The  requirement  to 
produce  and  store  radioactive  tritium  fuel  is  ehminated  and  the  net  inventory  of 
tritium  in  the  plasma  is  substantially  reduced.  This  reduces  the  associated  tritium 
handhng  problems  and  reduces  the  risk  of  serious  accidents  involving  tritium  releEise. 

The  disadvantages  of  the  D-^He  fuel  cycle  are  the  significant  physics  and  en¬ 
gineering  difficulties  which  must  be  overcome  to  make  D-®He  fusion  competitive. 
Relative  to  D-T  fusion,  D-^He  fusion  requires  higher  operating  temperature,  better 
confinement  of  the  plasma,  and  higher  plasma  beta.  Ignition  is  an  important  condi¬ 
tion  of  a  self-sustaining  fusion  reaction,  requiring  that  the  charged-particle  energy 
release  compensate  for  the  total  plasma  energy  loss  (particle  and  energy  transport 
and  radiation  losses).  Ignition  of  a  D-^He  plasma  in  terms  of  plasma  temperature 
T  (keV),  electron  density  and  plasma  energy  confinement  time  te  (s) 

is  optimized  at  TucTe  %  6  x  10”fceF  -  a  ’m"®  for  operating  temperature  T  55 
keV  [7].  For  D-T  ignition,  optimization  occurs  at  Tn^TE  «  3  x  lO^^keV  -  a 
for  operating  temperature  T  ss  15  keV,  substantially  lower  values  which  appear  to 
be  attainable  in  the  near  term.  Beta,  defined  as  /?  =  2fionT/B^  (i.e.,  the  ratio  of 
the  plasma  pressure  to  the  confining  magnetic  pressure),  is  an  important  parameter 
describing  tokamak  performsmce.  Beta  must  be  high  to  make  D-^He  fusion  practical 
since  the  fusion  power  density  at  equivalent  beta  values  is  much  lower  than  for  D-T 
[3].  Since  power  density  varies  eis  0^^  higher  betas  can  increase  it  to  be  comparable 
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with  D-T.  To  achieve  these  conditions  and  remain  economical,  however,  significant 
advances  in  fusion  reactor  physics  and  engineering  must  occur. 

An  additional  obstacle  to  competitive  D-®He  fusion  is  that  the  terrestriEd  sup¬ 
ply  of  the  helium  isotope  ^He  is  severely  Unoited.  While  the  hydrogen  isotope  deu¬ 
terium  is  found  in  plentiful  supply  in  nature,  terrestrial  ®He  resources  are  limited  to 
several  hundred  kilograms.  Assuming  a  consumption  rate  of  about  0.1  kg/MWe-3rr, 
these  resources  will  not  sustain  long-term  commercial  power  production.  Although 
it  may  be  possible  to  produce  ^He  as  a  b3rproduct  in  senoicatalyzed  deuterium  fusion 
reactors  [3],  in  which  aU  the  tritium  and  only  part  of  the  ^He  fuse  in  the  plasma, 
this  method  is  probably  not  cost  efiective.  Fortunately,  estimates  [8]  indicate  that 
the  regohth  of  the  lunar  surface  contains  about  a  million  tonnes  of  ^He  as  a  result 
of  deposition  by  the  solar  wind,  zmd  the  potential  quantity  of  ^He  in  Jupiter  and 
Saturn  is  about  10“  kilograms.  The  total  energy  cost  for  nouning  on  the  moon, 
degassing,  isotropic  separation  and  transportation  to  earth  is  estimated  [8]  to  be 
approximately  2.4  x  10®  MJ /kg  of  ®He.  The  thermal  energy  yield  from  D-®He  fusion 
is  about  6  x  10®  MJ/kg,  representing  ein  energy  payback  of  about  250  compared  with 
about  20  for  and  about  16  for  coal.  Although  an  economic  incentive  exists, 
msddng  lunar  supplies  of  ®He  practical  requires  a  new  and  expensive  industry,  de¬ 
pendent  upon  rehable  space  travel,  a  mining  operation  on  the  moon  and  numerous 
other  high-technology  advsmces.  Besides  the  physics  and  engineering  challenges  of 
reaching  actual  D-®He  ignition  2md  operation,  the  supply  of  ®He  itself  is  a  signifi¬ 
cant  obstacle  to  achieving  commercial  D-®He  fusion  power.  Therefore,  miTiiTig  the 
moon  and  the  outer  planets  for  reactsmt  fuels  such  as  ®He  would  require  further 
investigation  as  feasible  courses  of  action. 
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1.3  Tokamak  Fusion  Reactors 

The  tokamak  reactor  concept  gained  recognition  at  the  Novosibirsk  Conference 
in  1968  and  since  then  steady  progress  has  been  made  toward  its  development  into 
a  viable  fusion  core  for  a  thermonuclear  reactor.  Today  the  tokamak  is  the  most 
advanced  concept  for  magnetic  confinement  of  a  dense,  hot  plasma.  In  general,  it 
consists  of  a  toroidal,  axisymmetric  plasma  confined  by  the  combination  of  a  large 
toroidal  magnetic  field  (typically  5-10  T),  a  smaller  poloidal  magnetic  field  created  by 
a  toroidal  current  (typically  15-25  MA)  through  the  plasma  and  the  superposition  of 
equilibrium  and  shaping  fields  created  by  poloidal  field  coils  external  to  the  plasma. 
Figure  1.1  shows  a  schematic  of  the  basic  tokamak  configuration.  The  specific 
position  and  shape  of  the  plasma  cross  section  is  defined  by  the  magnetic  fields 
generated  by  the  poloidal  field  coils.  The  solenoidal  coil  in  the  center  hole  of  the 
tokamak  functions  as  a  transformer  to  induce  the  toroidal  current  in  the  plasma. 
There  is  a  technical  limitation  to  the  pulse  length,  determined  by  the  maximum 
magnetic  flux  that  can  be  provided  by  the  transformer  action,  making  the  toka.ma.lf 
2U1  inherently  pulsed  device  in  the  absence  of  some  external  current  drive  mechanism. 

The  Joint  European  Torus  (JET),  the  Tokamak  Fusion  Test  Reactor  (TFTR) 
at  Princeton,  DIII-D  in  San  Diego  and  JT-60  in  Japan  are  the  most  advanced 
operating  tokamaks  today.  The  International  Thermonuclear  Experimental  Reactor 
(ITER)  will  be  the  next  step  in  the  tokamzik  program  and  will  build  upon  the  opera¬ 
tional  experience  of  these  major  devices  and  other  tokamak  experiments  worldwide. 
These  world-class  tokamaks  are  large  machines.  JET,  for  example,  has  a  major 
plasma  radius  /2o=2.96m,  a  horizontal  minor  radius  a=1.25m,  and  a  vertical  minor 
radius  b=2.10m.  Figure  1.2  shows  the  JET  device  and  identifies  its  important  com¬ 
ponents,  t3rpic2d  of  most  current  tokamaks,  emd  Figure  1.3  shows  the  basic  geometry 
parameters  used  to  describe  tokamak  plasmas.  The  aspect  ratio  for  JET  is  A=2.37 
(an  important  parameter  defined  as  major  radius  Ro  divided  by  minor  radius  a)  and 
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pleisma  elongation  k=1.68  (defined  as  vertical  minor  radius  b  divided  by  horizontal 
minor  radius  a).  By  comparison,  the  aspect  ratio  for  TFTR  is  2.92,  for  DIII-D  is 
3.17  aud  for  JT-60  is  3.0. 

Results  from  tokamalc  experiments  over  the  past  decade  have  indicated  steady- 
state  operation  and  enhanced  confinement  of  tokamak  plasmas  may  be  achievable. 
Steady- state  operation  may  be  possible  by  noninductive  current  drive  provided  by 
high-energy  neutral  beam  injection  (NBI)  in  the  plasma  core,  radiofrequency  elec¬ 
tromagnetic  waves  in  the  outer  regions  of  the  plasma,  and  bootstrap  current  (i.e.  the 
current  generated  by  the  plasma  itsdf  -  see  A.4)  [9,  10].  The  advantages  of  steady- 
state  operation  include  reduction  of  thermal  fatigue,  mechanical  fatigue,  pla.sma 
disruption  and  reactor  down-time.  Enhanced  confinement  of  the  plasma  in  the  ”H- 
regime”  has  been  obtained  in  many  tokamaks  which  employ  a  magnetic  separatrix 
to  bound  the  plasma.  A  separatrix  configuration  with  a  plasma  divertor  and  high 
magnetic  shear  at  the  edge  is  thought  to  be  highly  advantageous  for  the  achieve¬ 
ment  of  high  confinement  H-mode,  which  offers  improvement  over  the  conventional 
Goldston  L-mode  confinement  scaling  by  a  factor  of  two  or  more  [11].  Achieving 
the  H-mode  requires  a  minimum  power  threshold  to  be  exceeded,  resulting  in  a 
bifurcation  with  instantaneous  changes  occurring  in  the  plasma  edge. 

For  H-mode  operation  in  a  separatrix  configuration  with  a  poloidal  divertor, 
the  effective  plasma  boundeiry  is  defined  by  the  separatrix  (also  referred  to  as  the 
plasma  X-point).  The  separatrix  surface  marks  the  boundary  between  the  closed 
magnetic  field  hnes  which  contain  the  plasma  and  the  open  field  hnes  which  leave 
the  vicinity  of  the  plasma  and  intercept  divertor  plates.  Tokamak  plasma  turbulence 
and  cross-field  transport  processes  cause  the  plasma  to  diffuse  across  the  internal 
flux  surfaces  and  eventually  across  the  separatrix  surface.  The  plasma  then  follows 
the  open  field  hnes  to  actively  cooled  divertor  plates,  which  absorb  the  kinetic 
energy  of  the  impacting  plasma  particles.  The  main  plasma-wall  interaction  occurs 
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at  these  plates,  where  the  plasma  is  neutralized  and  then  pumped  away  through 
nearby  ducts. 

Another  enhancement  of  plasma  confinement  is  thought  to  result  from  op¬ 
eration  in  a  second  stability  regime,  where  higher  beta  values  may  be  achievable. 
As  the  beta  values  increase,  the  plasma  becomes  susceptible  to  a  magnetohydro- 
d3nianiic  (MHD)  ballooning  instability.  Fluid  models  of  these  short  wavelength, 
localized  instabilities  predict  that  stability  is  sometimes  recovered  when  the  pres¬ 
sure  is  increased  above  that  corresponding  to  the  first  onset  of  instability.  It  is 
theoretically  possible  to  optimize  the  current  profile  in  such  a  way  that  ballooning 
instabilities  are  not  excited  at  any  beta  [11].  It  has  also  been  shown  numerically 
that  indentation  of  the  small-side  major  radius  may  result  in  stability  against  these 
MHD  modes  even  for  high  beta  values  and  may  provide  an  access  path  to  the  second 
stability  regime  [12]. 

Despite  the  enhancement  of  reactor  performance  resulting  from  steady-state, 
H-mode  operation  in  the  second  stability  regime,  the  attainment  of  beta  values  and 
plasma  confinement  good  enough  to  achieve  cost  effective  D-®He  fusion  may  be  con¬ 
strained  by  a  practical  lower  limit  on  the  eispect  ratio  for  the  basic  tokamak.  The 
conventional  wisdom  of  tokamadc  operation  and  prudent  engineering  suggests  the 
inclusion  of  a  solenoid  for  current-drive  start-up,  thick  nuclear  shields  inboard  to 
the  plasma  for  the  protection  of  superconducting  magnet  coils  and  insulators,  and  a 
separate  first  wall  and  vacuum  boundary.  These  factors  limit  the  major  radius  and 
aspect  ratio,  which  leads  to  modest  values  of  plasma  beta  (typically  around  5-10% 
for  aspect  ratios  of  around  3,  and  perhaps  greater  than  20%  in  the  2nd  stability 
regime)  and  require  high  on-axis  toroidal  magnetic  fields  for  plasma  confinement 
[13].  Plasma  beta  hmits  and  confinement  may  be  enhanced  with  high  plasma  elon¬ 
gation  [14],  but  to  achieve  plasma  elongation  of  even  1.6  with  mild  triangularity  for 
tokamaks  with  A  ~  3,  large  external  shaping  fields  (quadrupole  and  hexapole  fields) 
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are  required  [15].  The  high  capital  cost,  complexity,  and  maintenance  problems 
associated  with  a  large  fusion  power  core  and  superconducting  magnets  would  tend 
to  make  the  conventional  tokamak  reactor  design  an  unlikely  concept  for  practical 
D-^He  fusion. 

1.4  Spherical  Torus  Fusion  Reactor 

Reactor  studies  conducted  at  Oak  Ridge  National  Laboratory  in  1985  indicate 
performance  improvements  for  tokamaks  with  very  low  aspect  ratio  operating  at 
low  on-axis  magnetic  field.  Plasmas  with  a  very  low  aspect  ratio  {A  <  2)  display 
a  naturally  elongated  (k  >  2)  D-shaped  plasma  cross  section,  high  toroidal  beta 
(/3  >  20%),  and  other  desirable  features  [13,  16).  Since  the  magnetic  field  hnes  in 
a  very  low  aspect  ratio  plasma  have  a  high  pitch,  the  plasma  equilibrium  becomes 
highly  paramagnetic  and  essentially  force-free  with  a  large  plasma  current  density 
nearly  parallel  to  the  magnetic  field.  This  leads  to  a  strongly  enhanced  toroidal  field 
on-axis,  and,  combined  with  a  dominating  poloidal  field  at  the  outboard  region  of  the 
plasma,  creates  a  nearly  omnigeneous  region  which  is  nearly  free  of  locally  trapped 
particles,  contributing  to  the  kinetic  stability  of  the  plasma  and  reduced  neoclassical 
transport.  To  achieve  a  very  low  aspect  ratio,  inboard  components  such  as  the 
solenoid,  shield,  2uid  organic  insulator  are  eliminated  and  only  a  cooled,  normal 
conductor  that  carries  current  to  produce  the  toroidal  magnetic  field  is  retained. 
The  plasma  major  radius  is  significantly  reduced  and  the  resulting  plasma  looks 
much  like  a  sphere  with  a  small  hole  through  its  center,  hence  the  name  "spherical 
torus.”  Figure  1.4  shows  the  basic  design  configuration  of  a  sphericeil  torus  (ST) 
reactor.  The  features  of  a  spherical  torus  place  it  in  a  unique  physics  regime  which 
may  permit  compact  fusion  at  low  magnetic  field  and  modest  cost. 

Several  subsequent  studies  of  ST  reactors  have  investigated  ignition  assess¬ 
ments  and  engineering  feasibility  [17,  18,  19].  These  studies  indicate  ST  reactors 
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can  be  exceptionally  small  in  unit  size  and  highly  cost-effective,  but  axe  somewhat 
hmited  by  the  fact  that  no  experimental  data  base  for  very  low  aspect  tokamaks 
exist.  Assumptions  regarding  plasma  confinement  and  performance  are  extrapo¬ 
lated  from  the  current  physics  data  base,  accumulated  from  experiments  on  JET, 
TFTR,  DIII-D,  JT-60  and  other  research  reactors  (A  >  2.3),  which  may  or  may 
not  be  vahd  at  very  low  <ispect  ratio  (A  <  2).  Recently,  the  START  (Small  Tight 
Aspect  Ratio  Tokamak)  experiment  in  operation  at  Culham  Laboratory  has  formed 
hot  plasmas  at  aspect  ratios  around  1.3  which  exhibit  good  positional  and  MED 
stabihty,  and  exhibit  many  typical  tokamak  characteristics  [20,  21].  Oak  Ridge 
National  Laboratory  is  conducting  an  assessment  now,  and  preparing  to  produce 
a  conceptual  design,  for  a  hydrogen-fueled  steady-state  spherical  tokamak  (TST) 
intended  to  clarify  the  potential  of  the  TST  as  a  near-term,  low-cost  toroidal  device 
to  generate  critical  data  for  the  next  phase  of  magnetic  fusion  energy  development 
[22].  As  experimental  data  from  ST  reactors  become  available  to  verify  stability 
and  equilibrium  calculations  and  to  define  physics  scalings  appropriate  at  very  low 
Eispect  ratio,  more  credible  ST  reactor  studies  can  be  conducted. 

1.5  D-^He  ST  Fusion  Reactor 

Except  for  one  ST  reactor  study  investigating  the  ignition  and  bum  criteria 
for  a  D-®He  tokamak  [19],  eiU  ST  reactor  research  has  been  oriented  toward  D-T  fuel 
cycles.  The  current  ARIES-III  study  of  a  conceptual  D-®He  tokamak  reactor  assumes 
advanced  tokamak  parameters  and  advanced  physics  performance,  but  at  aspect 
ratio  around  3.0.  This  presented  an  opportunity  to  investigate  the  competitive 
potential  of  a  D-®He  ST  fusion  reactor,  which  was  the  basis  for  this  thesis.  E  compact 
D-®He  fusion  at  low  field  and  modest  cost  appears  viable,  then  the  possibility  may 
exist  for  competitive  fusion  power. 

Previous  D-T  ST  reactor  studies  conducted  at  Oak  Ridge  National  Laboratory 
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(ORNL)  used  a  computer  code  called  STORAC  (Spherical  TOrus  Reactor  Analysis 
Code),  run  on  the  Cray  computers  at  the  National  Energy  Research  Supercomputer 
Center.  This  code  is  a  ^Tiriant  of  TETRA  (Tokamak  Engineering  Test  Reactor  Anal¬ 
ysis),  an  earlier  reactor  analysis  code  developed  at  ORNL  which  calculates  tokamak 
performance,  cost,  and  configuration  as  a  function  of  plasma  and  engineering  pa¬ 
rameters  [23].  STORAC  assumes  a  0-dimensional,  profile  averaged  global  physics 
model  which  uses  conventional  tokamak  modeling  algorithms  modified  for  the  ST 
and  adopts  many  of  the  ITER  scalings  for  nominal  ST  mod^ng  (e.g.  energy  con¬ 
finement  laws,  etc.).  The  code  is  organized  in  a  standard  equation  solver  format, 
consisting  of  a  numerical  software  package  which  iterates  a  prescribed  set  of  vari¬ 
ables  in  order  to  satisfy  a  prescribed  set  of  constraints,  and  optionally  maximizes  or 
minimizes  a  given  figure  of  merit.  The  costs  accounting  used  in  STORAC  combines 
methods  used  in  TETRA  and  the  GENEROMAK  method  [24].  A  standeird  cost 
accounting  system  [25]  is  used,  which  follows  that  used  in  the  reporting  of  power 
reactor  costs,  updated  to  1990  dollars.  STORAC  is  benchmarked  with  ARIES-I  and 
ITER  computer  models  and  extrapolate  to  very  low  aspect  ratio  using  a  modified 
model. 

Using  STORAC  as  the  steirt  point,  the  research  for  this  theis  involved  esen- 
tially  five  steps.  They  were;  (1)  learning  the  architecture  of  the  existing  code  and 
how  to  run  the  code;  (2)  modifying  pints  of  the  code  to  convert  from  a  D-T  fuel  cycle 
to  a  D-®He  fuel  cycle;  (3)  debugging  the  modified  code  and  clarifying  assumptions; 
(4)  running  the  new  code  to  establish  credible  reference  cases;  and,  (5)  performing 
sensitivity  studies  for  several  key  parameters  of  the  reference  cases.  From  the  sensi¬ 
tivity  studies,  conclusions  regarding  the  potential  of  a  D-®He  ST  fusion  reactor  have 
been  made. 
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1.6  Thesis  Outline 

This  thesis  is  arranged  in  four  chapters  with  three  appendices.  Chapter  2 
reviews  the  modifications  made  to  STORAC  and  the  basic  assumptions  adopted  for 
this  study.  Chapter  3  discusses  the  reference  cases  and  the  results  of  the  sensitivity 
studies.  Chapter  4  discusses  the  economic,  safety  and  environmental  2issessment  of 
the  reference  cases  and  presents  the  conclusions  of  this  study.  Appendix  A  reviews 
the  spherical  torus  physics  model  used  in  STORAC.  (For  the  reader  unfamiliar 
with  fusion  reactor  physics,  a  brief  review  of  the  material  contained  in  Appendix  A 
would  be  beneficial  before  proceeding  to  the  next  chapter.)  Appendix  B  contains 
the  details  of  the  reference  cases  used  in  the  study.  Appendix  C  contains  a  copy  of 
all  the  output  plots  for  the  edge  safety  factor  scan  of  the  ST-2A  2nd  stability  regime 
reference  case,  which  is  representative  of  the  output  plots  produced  for  aU  the  other 
sensitivity  studies. 
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Coils  wound  around  torus  to 
produce  toroidal  magnetic  field 


Figure  1.1;  Schematic  of  a  tokamak 
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Figure  1.2:  The  JET  device.  (1)  Vacuum  vessel  (double  walled).  (2)  Material  limiter 
defining  the  outer  plasma  edge.  (3)  Poloidal  protective  shields  to  prevent  the  pleisma 
touching  the  vessel.  (4)  Toroidal  field  magnet  of  32  D-shaped  coils.  (5)  Mecheinical 
structure.  (6)  Outer  poloidal  field  coils.  (7)  Inner  poloidal  field  coils  (primary  or 
magnetizing  windings).  (8)  Iron  magnetic  circuit  (core  and  eight  return  sections). 
(9)  Water  and  electrical  connections  for  the  toroidal  field  coils.  (10)  Vertical  and 
radial  ports  in  the  vacuum  vessel  for  plasma  heating  systems  and  diagnostics. 


Figure  1.3:  Pl2isma  cross-sectional  shape  showing  major  radius  R,  horizonted  minor 
radius  a,  vertical  minor  radius  b,  triangularity  6  =  cf a,  and  indentation  d/2a.  The 
closed  lines  depict  magnetic  flux  surfaces. 


CHAPTER  2 

Spherical  Torus  Reactor  Analysis  Code  (STORAC)  and  Assumptions 

2.1  Introduction 

This  chapter  describes  the  computer  model  used  and  provides  justification  for 
the  assumptions  adopted  for  this  study.  Section  2.2  reviews  the  basic  architecture 
and  operation  of  the  original  code.  Section  2.3  describes  the  fundamental  modi¬ 
fications  made  to  the  code.  Section  2.4  reviews  the  assumptions  adopted,  which 
includes  justification  for  why  certain  constraints  and  variables  were  selected.  Meiny 
of  these  modifications  and  assumptions  were  developed  over  time,  based  on  several 
debugging  problems  and  preliminary  results  showing  fundamental  flaws  in  original 
assumptions.  The  modifications  and  assumptions  described  here  represent  those 
ultimately  adopted  for  this  study. 

2.2  Review  of  Original  Code 
2.2.1  Architecture  of  STORAC 

The  Spherical  Torus  Reactor  Analysis  Code  (STORAC)  consists  of  a  series 
of  modules,  each  describing  a  reactor  system  or  component,  controlled  by  a  driver 
or  optimizer  routine.  It  contains  a  numeric^tl  software  package  which  can  iterate  a 
prescribed  set  of  variables  to  satisfy  a  prescribed  set  of  constraints.  Tokamaik  reac¬ 
tors  are  complex  devices  consisting  of  many  nonhnear  interactions,  many  of  which 
are  described  in  Appendix  A.  STORAC  uses  the  generalized  nonlinear  program¬ 
ming  subroutine  VMCON  [26]  as  an  optimizer.  This  option  requires:  (1)  iterating 
more  variables  than  constraints  being  considered;  (2)  incorporating  upper  and  lower 
bounds  on  all  variables  being  iterated;  and,  (3)  specifying  a  prescribed  figure  of 
merit  to  be  minimized  or  maximized.  For  example,  using  cost  of  electricity  (COE) 
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as  the  figure  of  merit,  VMCON  is  able  to  find  the  device  configuration  providing 
minimum  COE  which  satisfies  the  prescribed  physics  and  engineering  constraints. 
Other  figure  of  merit  options  include  construction  cost,  aspect  ratio  and  plasma  ma¬ 
jor  radius.  Some  of  the  constraints  available  are  formulated  as  inequalities,  which 
provides  additional  flexibility  as  VMCON  iterates  to  find  a  self  consistent  set  of  de¬ 
vice  parameters.  An  option  to  use  a  nonlinear  equation  solver  without  optimization 
is  also  available,  which  is  useful  for  performing  benchmark  comparisons  using  fixed 
parameters. 

When  STORAC  runs  in  the  optimization  mode,  it  reads  in  an  input  file  identi¬ 
fying  which  constraints  are  to  be  used,  which  variables  are  to  be  used,  which  figure 
of  merit  is  to  be  used,  and  initial  values  for  the  various  device  parameters.  If  a 
parameter  is  not  identified  as  a  variable,  it  remains  fixed.  The  code  runs  through  its 
calculations  using  the  given  parameters.  The  VMCON  subroutine  then  calculates 
the  least  value  of  a  function  of  several  variables  subject  to  bnear  and  nonlinear  equal¬ 
ity  and  inequality  constraints.  Variables  are  then  incremented  upward  or  downward 
within  their  bounds  in  such  a  way  that  the  figure  of  merit  moves  closer  to  optimiza¬ 
tion.  The  code  continues  to  iterate  in  this  manner  until  all  prescribed  constraints  are 
met  to  within  a  given  error  tolerance,  such  as  10“'*.  When  optimization  is  complete, 
the  code  produces  an  output  file  listing  the  final  values  of  the  variables,  the  residujil 
error  of  the  constraints,  and  the  calculated  rea.ctor  parameters  which  include  system 
costs,  physics  parameter  values  and  engineering  parameter  values. 

2.2.2  D-T  ST  Model 

Before  addressing  the  modifications  made  to  STORAC  for  this  study,  it  is 
important  to  understand  the  D-T  ST  configuration  which  the  original  code  models. 
Basically,  the  D-T  fusion  reaction  at  optimum  operating  conditions  releases  about 
20%  of  its  17.6  MeV  reaction  energy  in  the  form  of  3.5  MeV  alpha  particles  and 


18 


about  80%  in  the  form  of  14.1  MeV  neutrons.  Other  fusion  reactions  are  negligible 
at  this  temperature  and  are  generaUy  ignored.  The  charged  alpha  particles  accu¬ 
mulate  as  ’’plasma  ash”  and  the  fzist  neutrons  escape  to  be  moderated  in  a  hthium 
blanket  surrounding  the  vacuum  vessel  and  first  wall  and  then  stopped  in  a  steel 
shield  surrounding  the  blanket.  The  lithium  blanket  is  important  for  breeding  tri¬ 
tium  for  later  use  as  fuel.  The  heat  energy  deposited  by  the  fast  neutrons  in  the 
blanket  and  shield  is  removed  by  a  coolant  fluid  and  converted  to  electricity  in  a 
basic  Rankine  cycle  emplo3dng  heat  exchangers,  turbine/generators,  condensers  and 
pumps.  Radiation  losses  from  the  fast  alpha  peirticles  are  small,  making  a  negligible 
contribution  to  the  overall  heat  energy  deposited  in  the  blanket /shield  structure  and 
available  for  conversion  to  electricity.  The  unshielded  center  conducting  post  is  also 
heated  by  the  fast  neutrons  and  must  be  cooled  to  keep  its  maximum  temperature 
below  a  threshold  for  material  melting  or  excessive  resistive  losses. 

The  D-T  ST  plasma  is  modelled  rdatively  simply  by  assuming  the  D-T  fuel 
ratio  Jind  alpha  ash  fraction  and  optimizing  the  plasma  electron  density  and  op¬ 
erating  temperature  for  a  spedfled  geometry.  (Many  other  parameters  must  also 
be  assumed  or  spedfled,  but  these  four  are  illustrative  of  how  the  code  operates.) 
The  neutron  fusion  power  resulting  from  the  D-T  reaction  is  converted  to  dectridty 
with  an  effidency  of  about  35-40%.  The  net  electric  power  produced  by  the  reactor 
equals  the  gross  electric  power  minus  the  resistive  losses  in  the  TF  magnet  coils 
(i.e,  the  center  conducting  post  and  TF  coil  return  legs),  plasma  heating  and  cur¬ 
rent  drive  power  requirements,  and  miscellaneouB  power  requirements  for  pumps, 
tritium  handling  systems,  etc.  Neutral  bejim  injection  is  generally  used  as  the  ex¬ 
ternal  current  drive  method,  assuming  a  deuterium/tritium  beam  interacting  with 
the  D-T  plasma  background.  For  a  given  energy  beam,  STORAC  calculates  the 
current  drive  effidency  as  well  as  the  number  of  beam-plasma  fusion  reactions  and 
the  beam  beta.  Based  on  the  plasma  physics  calculations  for  the  given  parameters, 
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STORAC  performs  engineering  calculations  to  determine  dimensions  of  the  reactor 
vessel,  CENTERPOST,  and  the  toroidal  field  (TF)  and  poloidal  field  (PF)  magnet 
coils  required  for  the  plasma  configuration.  Capital  costs  for  the  reactor  components 
are  then  calculated  and,  finally,  the  cost  of  electricity  is  determined. 

2.2.3  Constraint  Equations  and  Variables 

As  STORAC  iterates  to  optimize  a  prescribed  figure  of  merit,  certain  con¬ 
straints  are  necessary  to  insure  a  self  consistent  solution.  Constraint  equations 
which  must  be  satisfied  include:  the  plasma  beta  equation,  which  relates  beta  to 
the  plasma  pjirameters  (see  A.8);  the  global  plasma  power  balance,  which  relates 
power  sources  to  power  losses  (see  A.9);  the  beam  ion  density,  used  only  with  the 
NBI  option;  magnetic  field  at  the  cod,  which  insures  that  the  field  at  the  cod  is 
consistent  with  the  field  at  the  plasma;  the  radial  budd,  which  matches  the  inboard 
component  radial  budd  to  the  plasma  major  radius;  and  the  centerpost  temperature, 
which  insures  that  the  correct  resistivity  is  used.  Optional  inequality  constraints  re¬ 
lated  to  the  plasma  parameters  include:  the  density  limit  (see  A. 7);  the  Troyon  beta 
limit,  a  parameter  useful  in  first  stability  regime  calculations,  and  the  e/9p  limit,  a 
parameter  useful  in  second  stability  regime  calculations  (see  A.8);  the  volt  seconds 
limit,  which  is  used  for  pulsed  operation;  the  edge  q,  which  estabhshes  a  minimum 
edge  safety  factor  for  the  plasma  (see  A.3);  and  neutral  beam  energy,  where  beam 
decay  lengths  to  plasma  center  equals  a  prescribed  input  value.  Optional  inequality 
constraints  related  to  overall  device  operational  considerations  include:  the  neutron 
wall  load,  which  Hmits  the  £unount  of  radiation  damage  to  the  first  wall,  blanket 
and  shield;  the  divertor  heat  load,  which  limits  the  amount  of  erosion  damage  to 
divertor  plates  by  charged  particles;  the  fusion  power  limit;  the  net  electric  power 
hmit;  and  the  peak  centerpost  temperature. 

Numerous  physics  and  engineering  parameters  can  be  specified  as  variables, 
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with  prescribed  upper  and  lower  bounds.  Physics  parameters  include:  aspect  ratio 
A  =  Ro/s.]  plasma  major  radius  Rq',  edge  safety  factor  q^;  toroidal  magnetic  field  on- 
axis  Bq',  plasma  beta  0]  electron  density  11^;  plasma  temperature  Tg]  neutral  beam 
density  nbeam  ',  neutral  beam  energy  Ebeam',  confinement  H-factor;  and  several  inequal¬ 
ity  variables  for  use  with  the  physics  constraints.  Engineering  parameters  include: 
centerpost  current  density;  centerpost  average  temperature;  centerpost  coolant  ve¬ 
locity;  centerpost  coolant  fraction,  the  area  of  coolant  versus  cross  sectional  area 
of  the  centerpost;  TF  coil  thickness;  TF  coil  current  density;  and  several  inequality 
variables  for  use  with  the  operational  constraints. 

Initial  values  are  given  to  all  parameters  prescribed  as  variables  eis  wdl  as  all 
fixed  paremieters.  Fixed  parzimeters  include:  D-T  fuel  ratio;  alpha  ash  fraction; 
plasma  elongation  and  triangularity,  which  either  assume  a  fixed  scaling  dependent 
upon  inverse  aspect  ratio  (see  A.2)  or  axe  specified  explicitly;  Troyon  beta  scaling 
coefficient  and  maximum  e^p]  energy  confinement  time  scaling  (14  different  scaling 
options  available);  divertor  shape,  either  single  null  or  double  nuU;  divertor  model, 
either  conventional  or  expanded  geiseous  divertor;  plasma  scrape-off  layer  width  and 
length;  required  net  electric  power;  neutral  beam  ion  mass;  neutral  beam  width;  per¬ 
mitted  neutral  beam  e-decay  lengths  to  plasma  center;  bootstrap  current  fraction; 
angle  of  incidence  of  field  line  on  divertor  plate;  first  wall,  blanket  and  shield  thick¬ 
nesses;  centerpost  and  TF  coil  resistivities;  number  of  TF  coils  and  PF  coil  groups; 
coolant  inlet  temperature  and  thermal  conductivity;  thermal  to  electric  conversion 
efficiency;  and  unit  costs  and  other  cost  information  such  as  plant  capacity  factor, 
inflation  rate,  level  of  safety  assurance,  etc.  This  list  of  pareimeters  is  only  repre¬ 
sentative  of  the  many  parameters  required  to  run  the  code,  but  is  indicative  of  the 
complexity  of  a  tokamak  fusion  device.  It  is  important  to  note  that  initial  values 
used  to  start  the  code  should  be  reasonably  close  to  a  working  solution,  or  else 
the  physics  and  engineering  calculations  may  be  inconsistent  and  the  code  will  not 
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execute.  For  this  reason,  it  is  best  to  start  with  a  case  that  works  and  modify  the 
input  file  to  get  a  desired  output. 

2.2.4  Summary 

Although  this  review  of  the  original  code  is  brief,  it  provides  a  basic  description 
of  how  STORAC  runs,  what  it  is  trying  to  model,  and  some  of  the  constraints  and 
variables  available  for  modelling  a  reactor  configuration.  A  complete  description  of 
the  code  and  its  operation  is  beyond  the  scope  of  this  thesis,  which  is  not  meant  to 
be  a  STORAC  user  manual.  By  understanding  what  the  original  code  does,  one  will 
better  understand  what  the  modified  code  does  and  why  the  changes  were  made. 

2.3  STORAC  Modifications 
2.3.1  D-^He  ST  Model 

Identifying  the  general  features  of  the  D-®He  ST  configuration  was  necessary 
before  any  changes  could  be  made  to  STORAC.  How  would  the  D-®He  ST  be  dif¬ 
ferent  from  the  D-T  ST,  assuming  advanced  physics  parameters?  Conceptually, 
it  would  operate  at  low  field  in  steady-state,  with  high  beta  in  the  second  stability 
regime,  and  have  a  significant  plasma  current  dominated  by  the  bootstrap  effect  and 
supplemented  by  NBI  current  drive.  Tritium  fuel  would  not  be  required,  eliminating 
the  need  for  a  lithium  breeder  blanket  and  allowing  the  device  to  be  more  compact.  A 
direct  energy  convertor /expanded  gaseous  divertor  system  would  provide  conversion 
of  charged  particle  fusion  power  to  electricity,  offering  potentially  greater  efficiency 
than  the  standard  Rankine  cycle,  as  well  as  provide  plasma  ash  exhaust  and  im¬ 
purity  control.  Figure  2.1  shows  a  schematic  of  a  conceptual  double-null  poloidal 
divertor /direct  electrodynamic  energy  convertor  system  [6],  which  would  consist  of 
chambers  above  and  below  the  core  plasma  into  which  the  scrape-off  layer  plasma 
would  be  diverted,  collecting  ions  on  one  divertor  plate  and  electrons  on  the  other. 


22 


R(cm) 


Figure  2.1:  Schematic  of  poloidal  divertor/ direct  dectrodynamic  energy  conversion 
system. 

However,  in  the  D-^He  reactions  only  about  50%  of  the  total  charged*partide  energy 
rdease  is  subject  to  direct  conversion  due  to  radiation  and  high  paitide  transport 
to  the  walls  [7],  so  a  conventional  thermal  energy  conversion  system  would  still  be 
necessary  to  remove  a  significant  fraction  of  heat  energy  and  convert  it  to  dectridty. 

2.3.2  Code  Revisions 

The  first  significant  modification  required  to  transform  STORAC  from  a  D-T 
modd  to  a  D-^He  modd  was  to  incorporate  a  routine  to  solve  for  the  steady-state 
partide  balance  and  add  routines  for  additional  fusion  reaction  rates  (see  A.6).  A 
new  calculation  was  added  to  solve  for  the  deuterium  density  as  a  function  of 
partide  balance,  charge  balance  and  fud  ratio  nsfuD  and  then  determine  remaining 
partide  densities  as  a  function  of  n^.  At  the  optimum  D-^He  operating  temperature 
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in  the  50-80  keV  range,  the  D-T  reaction  rate  is  still  higher  than  for  D-*He  and 
significant  for  determining  neutron  fusion  power,  the  two  D-D  reaction  rates  are  low 
but  important  for  estabhshing  an  accurate  plasma  composition,  and  the  T-T  reaction 
rate  is  significantly  lower  but  accounts  for  some  removal  of  radioactive  tritium  from 
the  plasma.  Including  the  obvious  D-®He  reaction  rate,  these  other  four  reactions 
are  important  and  were  aU  considered  in  the  steady-state  particle  balance.  The 
original  code  simply  assumed  fractional  densities  of  deuterium,  tritium  and  alpha 
particles  relative  to  the  plasma  electron  density.  With  the  addition  of  ®He  and 
protons  to  the  plasma  background,  smd  with  five  different  fusion  reactions  occurring 
simultaneously,  no  simple  scaling  would  be  appropriate. 

The  new  particle  balance  required  new  variables  and  routines  to  be  defined 
and  old  ones  deleted  or  changed,  throughout  most  of  the  subroutines  in  the  code.  A 
new  variable  for  particle  confinement  time  Tp  required  a  new  constraint  equation  to 
be  satisfied  for  a  self  consistent  solution.  The  simple  D-T  scalings  for  charged  fusion 
power  and  neutron  fusion  power  (20%  and  80%  of  total  fusion  power,  respectively) 
had  to  be  generalized  with  new  calculations  for  charged  fusion  power,  neutron  fusion 
power  and  additional  fast  particle  losses  fusion  power.  The  D-T  scaling  for  fast 
alpha  beta  had  to  be  replaced  with  a  calculation  for  a  more  general  fast  particle 
beta,  accounting  for  all  fast  particles  produced  by  the  five  fusion  reactions  and  their 
different  energies.  The  calculations  for  fast  alpha  energy  transferred  to  electrons 
and  ions  were  also  replaced  with  new  calculations  for  generzd  fast  particle  energy 
transferred  to  the  electrons  and  ions.  The  NBI  calculations  had  to  be  modified 
to  account  for  a  D-®He  background,  essentially  by  replacing  tritium  entries  with 
®He  where  appropriate.  Since  the  operating  plasma  densities  for  D-®He  tend  to  be 
higher  than  D-T,  requiring  beam  energies  of  about  5  MeV  or  more  for  current  drive, 
the  NBI  calculations  had  to  be  modified  to  correct  a  fitting  factor  scaled  for  beam 
energies  of  about  1  MeV.  The  beam  fusion  calculations  also  had  to  be  modified  to 
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account  for  a  D-^He  background  and  new  beam  fusion  equations  added  to  account 
for  the  additional  reactions. 

The  next  significant  modification  to  STORAC  was  to  the  heating  power  and 
gross  electric  power  calculations.  Since  ®He  has  charge  Z=2,  the  electron  density  is 
a  larger  fraction  of  the  total  charged  particle  density  than  in  a  D-T  plasma.  This 
effect,  coupled  with  the  higher  operating  temperatures  for  D-^He,  means  radiation 
power  losses  play  a  larger  role  in  the  overall  power  balance  than  they  do  in  the  D-T 
pleisma.  Bremmstrahlung  and  synchrotron  radiation  losses  from  the  core  and  edge 
plasma  (see  A.9)  were  already  calculated,  but  neglected  from  the  equations  used  to 
determine  heat  energy  deposited  in  the  centerpost  and  shield  since  their  contribution 
was  negligible  compared  to  the  D-T  neutron  energy.  These  equations  were  modified 
to  include  both  neutron  and  radiation  energy  deposited  in  the  centerpost  zind  shield. 
The  inclusion  of  a  direct  energy  convertor  required  an  additional  calculation  to 
convert  charged  partide  power  to  the  divertor  into  dectridty  vrith  a  prescribed 
DEC  effidency.  The  remaining  power  to  the  divertor  would  then  appear  eis  heat 
energy  to  the  divertor  plates,  to  be  induded  in  the  heat  energy  sum.  The  total 
heat  energy  deposited  to  the  shield  and  divertor  plates  would  then  be  converted  to 
dectridty  with  a  prescribed  thermal  effidency.  Thus,  the  gross  dectric  power  would 
be  a  sum  of  the  dectridty  generated  by  the  DEC  and  the  turbine/generators. 

2.3.3  Summary 

Adding  the  steady-state  partide  balance  routine  to  the  code,  with  its  assod- 
ated  modifications  to  other  variables  Jind  routines,  and  modifying  the  gross  dectric 
power  calculations  to  account  for  direct  energy  conversion  represented  the  significant 
modifications  to  STORAC.  Most  of  these  changes  were  obvious  from  the  start,  while 
others  became  evident  once  the  code  was  running  and  produdng  output  of  dubious 
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quality.  To  list  here  each  individual  code  change  made  to  STORAC  would  be  un¬ 
necessarily  tedious.  (The  individual  code  changes  are  documented  in  the  STORAC 
files.)  What  is  important  to  understand  here  are  the  fundamental  modifications  to 
STORAC  that  were  necessary  for  modelling  the  new  D-^He  ST  configuration. 

2.4  Assumptions  Made  for  this  System  Study 
2.4.1  Introduction 

With  the  code  modifications  to  STORAC  in  place,  the  task  of  establishing 
reference  cases  for  further  study  required  a  number  of  assumptions  oriented  toward 
specific  objectives.  Since  the  stated  purpose  of  this  study  was  to  investigate  the 
potential  of  the  steady-state  D-^He  ST  for  commercial  electric  power  production,  the 
reference  cases  should  seek  to  minimize  the  cost  of  electricity  for  a  typical  1000  MWe 
fusion  reactor.  One  case  should  assume  operation  in  the  first  stability  regime  and  one 
case  should  aggressively  assume  operation  in  the  theoretical  second  stability  regime. 
Both  cases  should  eissume  advanced,  though  reasonable,  physics  smd  engineering 
parameters.  Once  these  two  reference  cases  are  established,  sensitivity  studies  can 
investigate  the  effect  of  varying  the  values  of  important  parameters. 

The  reference  cases  for  this  study  both  used  13  constraint  equations  and  21 
variables  to  arrive  at  a  self  consistent  successful  optimization.  The  following  con¬ 
straints  were  used  to  develop  these  cases:  plasma  beta;  global  plasma  power  balance; 
neutral  beam  ion  density;  field  at  coil;  centerpost  temperature;  particle  confinement 
time;  net  electric  power  limit;  and  neutral  beam  energy.  The  following  inequality 
constraints  were  used  to  develop  the  reference  cases:  density  limit;  Troyon  beta 
limit;  ej3p  limit;  and  edge  safety  factor  q.  The  following  variables  axe  used:  aspect 
ratio,  A=i2o/a;  toroidal  magnetic  field,  Bq\  plasma  major  radius,  Roi  average  elec¬ 
tron  temperature,  (Te);  plasma  beta,  /?;  average  electron  density,  (rie);  neutral  beam 
density  fraction,  (rie);  confinement  H-factor,  /j,;  current  drive  power  used  for 
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additional  heating;  neutral  beam  energy,  Ey^arni  centerpost  radius;  centerpost  cur¬ 
rent  density,  J^p]  centerpost  average  temperature;  centerpost  coolant  velocity;  aver¬ 
age  radius  of  centerpost  coolant  channel;  centerpost  coolant  fraction;  and  TF  coil 
current  density.  Four  other  variables  were  used  for  the  inequahty  constraints.  The 
figure  of  merit  to  be  minimized  was  cost  of  electricity.  The  residual  error  tolerance 
for  the  VMCON  optimizer  was  10“^. 

2.4.2  Discussion  of  Variables 

The  basic  geometry  of  the  plasma  and  the  reactor  required  to  maintain  it 
is  dependent  upon  the  aspect  ratio  and  plasma  major  radius.  The  aspect  ratio 
was  allowed  to  vary  between  a  lower  limit  of  1.2  [18]  and  an  upper  Hmit  of  2.0, 
representing  the  full  operational  regime  of  ST  plasmas.  The  plasma  major  radius 
was  allowed  to  vary  over  a  broad  range  of  1-4  meters.  These  Hmits  would  allow  for 
optimization  at  whatever  plasma  size  is  required  for  a  1000  MWe  reactor,  although 
both  aspect  ratio  and  major  radius  tend  to  be  as  low  as  possible  for  cost  effective 
performance. 

The  plasma  performance  is  dependent  upon  operating  temperature  and  den¬ 
sity,  as  well  as  the  beta  value  required  for  stable  confinement  and  the  confinement 
H-factor  required  for  power  balance.  Since  D-®He  fusion  peaks  at  temperatures  of 
about  55  keV,  the  electron  temperature  was  given  a  reasonable  range  of  40-80  keV. 
Density  was  assumed  to  vary  between  10*°  —  10*^m“*,  with  the  Borass  density  hmit 
(see  A. 6)  expected  to  be  within  this  range  and  the  density  to  be  at  or  sUghtly  below 
this  hmit  for  maximum  fusion  power  density.  Since  strong  paramagnetism  intro¬ 
duces  an  important  uncertainty  in  the  apphcation  of  the  Troyon  beta  hmit  [16], 
the  plasma  beta  was  given  a  broad  range  of  0.20-2.0,  with  the  upper  hmit  set  high 
enough  to  allow  for  possible  beta  values  greater  than  1.0  for  operation  in  the  second 
stabihty  regime  [27].  The  confinement  H-factor  was  also  given  a  broad  range  of  1-10 
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to  eiUow  for  whatever  H-mode  enhancement  was  required  for  global  plasma  power 
balance  (see  A.9).  All  of  these  bounds  were  broad  enough  so  as  to  not  constrain  an 
optimization  unnecessarily. 

The  magnetic  confinement  of  the  plcisma  is  dependent  on  the  magnet  coil  con¬ 
figuration  surrounding  the  reactor  vessel.  The  toroidal  magnetic  field  on-axis  was 
allowed  a  range  of  1-4  Tesla,  with  the  field  strength  expected  to  fall  in  the  2-3  T 
range  for  the  assumed  reactor  configuration  [27].  To  produce  the  required  mag¬ 
netic  field  on- axis,  strong  electric  currents  must  pass  through  the  centerpost  2md 
TF  coil  legs.  The  allowable  current  density  through  the  centerpost  is  limited  by 
temperature  and  material  stress  constraints  (however,  only  the  temperature  con¬ 
straint  is  used  in  the  current  model).  Since  the  maximum  allowable  current  density 
determines  the  minimum  radius  of  the  centerpost,  to  which  the  aspect  ratio  is  very 
sensitive,  the  centerpost  parameters  are  critical.  Assunoing  an  imshidded  center- 
post  could  be  made  of  a  single  dispersion  strengthened  copper  conductor  cooled  by 
high-velocity  pressurized  water,  the  maximum  current  density  is  estimated  to  be 
about  75  MA/m*  without  causing  the  peak  copper  temperature  to  exceed  250®C 
[18].  Thus,  the  centerpost  current  density  w£is  given  an  upper  bound  of  75  MA/m*, 
the  average  temperature  was  given  the  upper  bound  of  125  °C,  Md  the  centerpost 
radius  was  allowed  a  range  of  0.1- 1.0  m.  The  coolant  velocity  of  high-speed  pressur¬ 
ized  water  was  assumed  to  have  an  upper  bound  of  20  m/s  [17],  which  neglects  the 
corrosive  effects  of  high  velocity  coolant.  The  average  radius  of  the  coolant  channels 
was  allowed  a  range  of  0.001-0.010  m  and  the  coolant  fraction  of  the  centerpost  was 
allowed  a  range  of  0. 1-0.4,  where  the  number  of  coolant  channels  is  calculated  to 
equal  the  coolant  area  divided  by  the  average  coolant  channel  area.  For  the  con¬ 
ventional  copper  TF  coil  legs,  which  do  not  involve  such  critical  parameters,  the 
current  density  was  allowed  a  range  of  1-10  MA/m*.  The  dimensions  of  the  TF  coils 
are  then  determined  by  dividing  the  required  current  by  the  current  density  to  get 
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the  cross  sectional  area,  which  is  then  multiphed  by  the  required  TF  leg  length  to 
give  the  volume  of  the  TF  coil  legs.  STORAC  optimizes  to  find  the  best  bEilance 
between  the  cost  of  the  centerpost  and  TF  coil  magnet  assembly  (volume  x  unit 
costs)  and  resistive  power  losses,  based  on  all  prescribed  parameters. 

The  noninductive  current  drive  provided  by  NBI  requires  high  energy  neutral 
beams  to  generate  an  ion  current  equal  to  the  injected  beam  current  multiplied  by 
the  average  number  of  times  the  fast  ions  circulate  the  torus  before  being  thermahzed 
[11].  Counteracting  this  current  is  a  reverse  electron  current  due  to  friction  with 
the  bulk  electrons,  which  forms  a  component  of  the  plasma  current.  This  generated 
current  is  dependent  upon  the  density  and  energy  of  the  neutral  beam.  STORAC 
determines  the  beam  density  and  energy  required  for  current  drive.  The  beam 
density  as  a  fraction  of  the  electron  density  was  allowed  a  range  of  10  —  10  and 

the  beam  energy  was  allowed  a  range  of  1-10®  keV  to  allow  full  flexibility  in  meeting 
current  drive  requirements.  The  neutral  beam  power  used  for  additional  plasma 
heating  (if  reqmred)  was  allowed  a  range  of  1-100  MW. 

The  remaining  four  variables  were  used  for  the  density  Hmit,  Troyon  beta  Hmit, 
e/?p  limit,  and  minimum  edge  q  inequality  constraint  equations,  with  their  respective 
’f- values’  allowed  a  range  of  0. 0-1.0.  These  variables  would  allow  the  average  electron 
density,  plasma  beta  and  e/Sp  values  to  be  at  or  below  their  upper  limits  and  for  the 
edge  q  value  to  be  at  or  above  its  lower  limit.  This  would  allow  for  considerable 
flexibihty  during  optimization  since  each  and  every  constraint  would  not  have  to  be 
at  an  extreme. 

2.4.3  Discussion  of  Fixed  Parameters 

This  intent  of  this  section  is  to  identify  most  of  the  fixed  parameters  used  to 
estabhsh  the  reference  cases,  with  some  justification  for  why  they  were  chosen.  Many 
of  the  parameters  used  as  defaults  in  the  physics  model  are  described  in  Appendix 
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A  and  will  be  summaiized  here.  Other  physics  and  engineering  parameters  critical 
for  this  study  will  also  be  described  in  this  section. 

Plzisma  shape  and  geometry  (see  A.2)  assumed  plasma  elongation  and  trian¬ 
gularity  to  scale  as  a  function  of  inverse  aspect  ratio.  This  study  also  eissumed  a 
plasma  scrape-oif  layer  thickness  of  0.04  m,  a  first  wall  thickness  of  0.02  m,  and  a 
shield  thickness  of  0.50  m.  Since  no  blanket  is  required,  the  bleinket  thickness  was 
set  equal  to  zero.  The  vertical  gap  between  the  shield  and  TF  coils  was  assumed 
to  be  0.163  m.  A  double-null,  expanded  gaseous  divertor  concept  wa^  assumed,  in 
which  a  simple  divertor  heat  load  calculation  is  employed  (i.e.,  the  divertor  power 
is  assumed  to  be  uniformly  distributed  about  the  divertor  area  by  the  action  of 
copious  amounts  of  neutral  particle  interactions).  The  increase  in  the  vertical  build 
to  accommodate  this  divertor  concept  was  then  accounted  for  in  the  vertical  build. 

Other  shape  and  geometry  pareimeters  were  related  to  the  magnet  parameters. 
The  conventional  copper  TF  coil  legs  were  assumed  to  have  a  resistivity  of  2.5x10® 
Ohm-m  and  the  strengthened  copper  centerpost  was  assumed  to  have  a  resistivity 
of  3.0x10®  Ohm-m.  The  reactor  is  assumed  to  be  built  on  top  of  or  adjacent  to 
its  power  supplies,  with  a  bus  length  of  10  m  for  each  leg  and  the  bussing  current 
density  aissumed  to  be  1.25x10®  A/m*.  Two  superconducting  PF  coil  groups  are 
assumed,  with  the  PF  coils  assumed  to  be  external  to  the  TF  coils. 

The  assumed  plasma  current  scaling  (see  A.3)  is  based  on  a  double-null  divertor 
scaling.  The  minimum  edge  q  was  scaled  as  a  function  of  inverse  aapect  ratio.  The 
maximum  bootstrap  current  fraction  (see  A.4)  was  assumed  to  be  0.9.  NBI  provides 
the  remaining  noninductive  current  drive,  which  provides  the  seed  current  for  the 
bootstrap  effect  and  profile  control  for  the  plasma  current.  The  beam  width  wsis 
assumed  to  be  0.31  m  and  the  neutral  beam  wall  plug  to  injector  efficiency  was 
assumed  to  be  0.5.  The  permitted  neutral  beam  e-decay  length  to  plasma  center 
ratio  was  assumed  to  be  3.0  (i.e.,  nj,eam(0)  =  e~®nbeom)- 
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Other  miscellaneous  physics  parameters  are  identified  in  Appendix  A.  The 
density  and  temperature  profile  exponents  (see  A.5)  were  assumed  to  be  0.5  and 
1.0,  respectively.  The  ion  and  electron  temperatures  were  assumed  to  be  equal.  The 
particle  confinement  time  to  energy  confinement  time  ratio  (see  A.6)  was  assumed 
to  be  3.0.  The  effective  helium-3  to  deuterium  fuel  ratio  (see  A.6)  was  assumed  to 
be  1.0.  Tritium  recycling  was  assumed.  The  trapping  fractions  of  the  fast  charged 
particles  resulting  from  fusion  reactions  (see  A.6)  were  assumed  to  be  1.0,  with  the 
exception  that  0.9  was  assumed  for  the  14.7  MeV  protons.  The  1990  ITER  H-mode 
power  scaling  (see  A. 9)  was  assumed  for  the  energy  confinement  time.  The  reactor 
wall  reflectivity  (see  A.9)  was  assumed  to  be  0.9,  representative  of  a  hot  stainless 
steel  vacuum  chamber  at  400 °C  [28]. 

Among  the  remaining  engineering  parameters,  the  average  thermal  energy 
conversion  efficiency  was  assumed  to  be  40%  (t3rpical  of  most  power  plants)  and  the 
direct  electrodynamic  conversion  effidency  was  assumed  to  be  50%  [6].  The  levd 
of  safety  assurance  (LSA)  was  assumed  to  be  2,  which  assumes  large-scale  peissive 
protection  in  the  event  of  an  acddent  [1].  The  plant  operating  life  £md  capadty  factor 
were  assumed  to  be  30  years  and  0.75,  respectively.  The  cost  accounting  remained 
the  same  as  established  in  the  origineil  code  [24,  25],  except  for  two  modifications 
[27].  The  unit  cost  for  the  divertor  system  was  assumed  to  be  $1  million/m*  to 
account  for  the  additional  complexity  of  the  direct  energy  conversion  system  and 
the  unit  cost  for  the  NBI  system  was  assumed  to  be  $5 /watt  to  account  for  high 
energy  beam  systems  (greater  than  1  MeV). 

The  reference  cases  for  the  first  2md  second  stabihty  regimes  (see  A. 8)  assumed 
the  same  parameters  except  three.  For  the  first  stabihty  regime  reference  case,  the 
Troyon  coeffident  for  the  beta  hmit  was  eissumed  to  be  3.3  for  steady-state  operation, 
the  efip  hmit  was  assumed  to  be  0.6,  and  the  central  safety  factor  q(0)  was  assumed 
to  be  1.0.  For  the  second  stabihty  regime  reference  case,  the  Troyon  beta  hmit  was 
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not  used,  the  ej0p  limit  was  assumed  to  be  1.0,  and  q(0)  was  assumed  to  be  2.0  [27]. 
All  other  par2uneters  discussed  previously  remained  the  same. 

2.4.4  Summary 

Having  decided  upon  which  constraint  equations  and  variables  to  use,  with 
appropriate  bounds  and  initial  values,  the  task  of  estabhshing  reference  cases  was 
reduced  to  running  the  code  to  get  a  successful  optimization  eind  checking  that  the 
output  variable  values,  when  re-entered  as  initial  values  in  the  input  file,  converged 
to  a  self-consistent  solution.  Based  upon  all  assumptions  made,  there  should  be  only 
one  solution  for  each  reference  case.  Any  change  to  the  input  file,  such  as  assuming 
a  thermal  conversion  efficiency  of  only  35%  rather  than  40%,  would  specify  a  dif¬ 
ferent  reactor  and  result  in  a  difierent  solution.  For  the  most  part,  the  assumptions 
defined  here  were  agreed  upon  [27]  throughout  the  period  of  code  modification  and 
debugging. 


CHAPTER  3 

Results  of  D-^He  ST  Reactor  Study 


3.1  Introduction 

This  chapter  discusses  the  reference  cases  which  resulted  from  the  assumptions 
described  in  the  previous  chapter  and  the  results  of  several  parametric  investigations. 
Section  3.2  discusses  the  reference  cases,  both  the  initial  results  and  the  final  ref¬ 
erence  cases  selected  for  further  study.  Section  3.3  investigates  the  sensitivity  of 
the  2nd  stability  reference  case  to  efip,  which  was  given  a  range  of  0.5-1.5.  This 
is  primarily  to  demonstrate  a  linkage  between  the  1st  and  2nd  stabihty  regimes. 
With  a  linkage  between  the  two  reference  cases  established,  further  investigations 
can  focus  on  just  one  reference  case,  that  being  the  more  optimistic  2nd  stability 
case.  The  remaining  sections  investigate  the  sensitivity  of  the  2nd  stability  reference 
case  to  the  following  nine  parameters;  (1)  temperature  and  density  profile  factors; 
(2)  ratio  of  particle  confinement  time  to  energy  confinement  time;  (3)  fuel  ratio; 
(4)  confinement  H-factor;  (5)  bootstrap  current  fraction;  (6)  aspect  ratio;  (7)  cen- 
terpost  current  density;  (8)  thermal  conversion  efficiency;  and,  (9)  DEC  conversion 
efficiency.  Reactor  performance  based  upon  the  varying  parameter  values  is  then 
quantitatively  evaluated  relative  to  the  2nd  stability  reference  case. 

3.2  1st  and  2nd  Stability  Regime  Reference  Cases 
3.2.1  Initial  Results 

Based  upon  the  basic  assumptions  adopted  for  this  study,  STORAC  produced 
successful  optimizations  for  a  1st  stability  regime  ST  reactor  {efip  =  0.6)  and  a  2nd 
stability  regime  ST  reactor  (e/3p  =  1.0).  These  two  cases  have  been  designated  ST- 
lA  and  ST-2A,  respectively.  Table  3.1  lists  the  key  parameters  of  these  two  cases. 
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This  table  also  lists  the  key  parameters  for  ARIES-III  [29],  a  conceptual  D-®He 
2nd  stability  regime  tokamaJc  reactor,  as  a  benchmark  comparison.  Although  the 
ARIES-III  results  are  about  as  speculative  as  the  results  of  this  study,  considerable 
expertise  has  been  committed  to  the  ARIES  program  and  the  ARIES-III  results 
provide  a  useful  benchmark  for  evaluating  the  relative  potential  of  the  D-®He  ST 
reactor. 

The  bottom  Hne  comparison  of  these  results  indicate  that  the  estimated  unit 
cost  of  electricity  (COE)  in  miU/kWh  (mills  per  kilowatt-hour;  1000  Tuills  =  $1)  is 
116.5  for  ST-lA,  88.5  for  ST-2A,  and  74.3  for  ARIES-III.  Costs  are  in  1990  doll2urs 
for  all  three  cases.  Comparable  COE  estimates  [30]  for  an  advanced  pressurized 
water  reactor  (APWR)  and  a  soIm  photovoltaic  electric  plant  are  40-50  mill/kWh 
and  200-300  miU/kWh,  respectively.  Thus,  the  D-®He  fusion  reactors  do  not  appear 
to  be  as  economically  competitive  as  advanced  fission  breeder  reactors  (assuming 
COE  comparable  with  APWR),  but  they  do  appear  to  fall  within  the  broad  COE 
range  expected  for  future  energy  systems. 

Comparison  of  ST-IA,  ST-2A  eind  ARIES-III  results  indicate  significantly 
lower  toroidal  magnetic  fields  on-axis  for  the  ST  reactors,  as  expected.  However, 
the  higher  peak  magnetic  fields  required  at  the  TF  coil  for  the  ST  reactors  are 
characteristc  of  very  low  aspect  ratio  geometry.  This  may  be  seen  in  the  following 
equation,  where  the  magnetic  field  at  the  plasma  axis  is 


Bo  —  Be 


Rq  ~  tto  —  * 
Rc 


where  s  is  the  distance  between  the  plasma  and  magnet  coil  necessary  to  house  the 
first  wall  and  shield.  The  higher  fields  required  at  the  TF  coil  lead  to  higher  resistive 
power  losses  in  the  copper  magnets,  resulting  in  Isirger  recirculating  power  fractions 
and  higher  COE.  Ideally,  the  reference  cases  for  the  ST  should  have  recirculating 
power  fractions  comparable  with  ARIES-III,  where  20-30%  represents  a  reasonable 
range.  Power  losses  much  greater  than  300-400  MWe  for  a  1000  MWe  commercial 
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power  plant  simply  are  not  attractive  to  the  dectric  utilities. 

A  few  other  points  are  worth  noting  from  Table  3.1.  First,  the  neutron  wall 
load  is  much  lower  for  ARIES-III  than  for  the  ST  reactors,  since  the  first  wall  surface 
area  for  ARIES-III  is  larger  (1394  m*  versus  less  than  400  m^  for  the  ST  reactors). 
Second,  the  total  magnet  costs  are  also  lower  for  ARIES-IU  than  for  the  ST  reactors, 
since  the  ST  reactors  use  normal  copper  magnets  and  more  than  60%  of  the  toted 
magnet  costs  are  for  the  PF  magnet  assemblies,  which  are  necessary  for  the  stronger 
vertical  fields  required  for  ST  confinement.  Third,  the  NBI  beam  energy  for  ARIES- 
III  appeeirs  to  be  very  large,  but  the  models  used  in  both  ARIES  and  STORAC  are 
very  speculative.  The  beam  energies  and  power  requirements  are  useful  for  cost 
estimates,  but  further  conclusions  should  be  reserved  until  better  models  become 
available.  Lastly,  the  Lawson  parameter  ute  (which  defines  the  breakeven  condition 
for  a  fusion  plasma)  for  both  ST-IA  and  ST-2A  is  about  1.8  x  10*^s/to®,  which 
appears  to  be  well  within  the  ignition  xepme.  Previous  studies  [19]  indicate  D-*He 
ignition  in  a  spherical  torus  could  be  achieved  with  wall  reflectivities  below  80%  and 
at  ute  <  10*^s/m®,  without  strong  plasma  shaping  or  e0p  >  0.6. 

3.2.2  Variation  of  Edge  Safety  Factor 

The  puameter  with  perhaps  the  most  significant  effect  on  reactor  performance 
was  the  edge  safety  factor,  defined  as  =  aBo/RoBp.  The  edge  sedety  factor  seeding 
used  for  the  q-limit  constraint  equation  is  aspect  ratio  dependent  (see  A.3),  which 
may  or  may  not  be  valid  for  very  low  eispect  ratio.  Theoretically,  a  high-q  reactor 
seems  adveintageous  with  regard  to  decreasing  the  statistical  probability  of  plasma 
disruptions  [9].  Disruptions  due  to  MED  instabilities  can  lead  to  the  sudden  loss  of 
plasma  confinement,  an  operational  transient  which  could  residt  in  an  accident.  In 
addition  to  increased  risk  of  disruptions,  plasma  confinement  and  stabihty  also  tend 
to  deteriorate  with  low-q  operation.  Although  ~  3  appears  to  be  the  lower  limit 
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for  safe  operation,  the  practical  lower  limit  is  probably  higher  and  only  detailed 
MHD  calculations  for  a  specific  reactor  configuration  could  determine  what  the 
actual  q-limit  should  be. 

Reducing  the  edge  safety  factor  reduces  the  on-axis  toroidal  magnetic  field  and 
the  required  peak  field  at  the  TF  coil,  leading  to  lower  resistive  power  losses  and 
lower  COE.  Figures  3.1  through  3.4  show  various  effects  on  reactor  performance  as 
is  reduced  from  7.68  to  3.5  for  ST-IA  and  ST-2A,  respectively.  As  seen  in  the 
figures,  there  is  a  linear  reduction  in  the  on-axis  magnetic  field  Bo  with  lower 
values,  leading  to  sm2iller  major  radius  values  and  much  smaller  TF  coil  resistive 
power  losses.  The  required  injection  power  also  shows  a  linear  decrease  with  lower  q 
values.  The  lower  recirculating  power  fractions  associated  with  lower  values  mean 
less  fusion  power  is  required  to  meet  the  net  electric  power  constraint,  resulting  in 
more  compact  reactors  and  lower  costs. 

3.2.3  Reference  Cases  Used  for  the  Study 

Based  on  the  results  of  low-q  effects  on  reactor  performance,  new  reference 
cases  for  the  study  were  estabhshed  which  attempt  to  strike  a  balance  between  the 
physics  advantages  of  high-q  operation  (i.e.,  reduced  risk  of  disruptions)  and  the 
engineering  advantages  of  low-q  operation  (i.e.,  reduced  resistive  power  losses  and 
lower  COE).  An  edge  safety  factor  =  5  was  selected  as  a  reasonable  value  on  which 
to  base  further  parametric  investigations  [27].  Since  high  0  hmits  and  good  energy 
confinement  are  strongly  correlated  to  large  plasma  current,  which  leads  to  operation 
at  low  q,  values  of  q^  <  5  may  be  attainable.  Table  3.2  lists  the  key  puameters  of 
these  two  reference  cases,  designated  ST-IB  and  ST-2B,  respectively.  Appendix  B 
contedns  a  complete  output  listing  of  all  pareimeters  for  these  two  reference  cases. 

The  cost  of  electricity  for  both  ST-IB  and  ST-2B  (88.9  and  73.7  mill/kWh, 
respectively)  are  comparable  with  ARIES-III,  due  to  the  reductions  in  reactor  size. 
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magnetic  Adds  and  recirculating  power  fractions  associated  with  a  lower  q-value. 
The  use  of  superconducting  magnets  and  lower  PF  Add  requirements  in  ARIES- 
III  tends  to  offset  the  cost  of  high  Add  operation,  which  the  ST  reactors  must 
compensate  for  with  better  plasma  performance  at  very  low  aspect  ratio.  The  higher 
beta  values  permit  lower  Add  requirements  and  higher  operating  temperatures  and 
densities,  which  lead  to  fusion  power  densities  high  enough  to  make  ST  reactors 
competitive  with  advanced  tokamaks.  Both  reference  cases  have  a  major  radius 
of  about  2  m  and  an  on-ajds  magnetic  Add  of  less  than  2  T.  Despite  the  rdative 
’compactness’  of  these  two  cases  {Rq  2  m),  these  are  not  small  machines;  the  radial 
build  for  both  cases  is  about  5  m  and  the  vertical  build  is  more  thcin  8m-  roughly 
30’  in  diameter  and  50’  high.  The  total  capital  cost  for  both  cases  is  estimated 
at  about  $5  billion,  which  would  be  a  signiAcant  investment  for  any  dectric  power 
utility. 

As  indicated  in  Chapter  2,  all  assumptions  used  to  devdop  these  reference 
cases  are  based  upon  reasonable  extrapolations  of  physics  knowledge  and  engineering 
technology.  SpedAc  details  will  have  to  be  resolved,  regardless  of  siny  assumptions 
adopted  for  this  study.  Although  the  spedAc  cost  accounting  used  in  STORAC 
may  be  different  from  that  used  in  the  ARIES-III  study  (unit  costs  assumed  for 
ARIES  appear  to  be  lower  than  those  assumed  for  this  study),  a  useful  benchmark 
comparison  can  be  made  between  the  ST  reference  cases  and  ARIES-III.  Therefore, 
these  reference  cases  describe  reeisonably  well  the  antidpated  characteristics  of  future 
D-^He  ST  reactors,  upon  which  useful  parsimetric  investigations  may  be  conducted. 

3.3  Variation  of 

Although  it  has  been  observed  that  gross  plasma  conAnement  in  1st  stabil¬ 
ity  regime  tokamaks  deteriorates  for  e0p  >  0.6,  practiced  limits  on  e/9p  values  for 
operation  in  the  2nd  stability  regime  are  stiU  undeAned.  Values  approaching  unity 
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are  believed  to  be  reasonable,  but  ARIES-III  assumes  e^p  =  1.803,  suggesting  that 
higher  values  may  be  attainable.  Figures  3.5  and  3.6  show  the  effect  on  the  ST-2B 
reactor  performance  for  eySp  values  over  the  range  0.5-1.5.  As  seen  in  these  fig¬ 
ures,  reactor  performance  generally  improves  with  higher  efSp.  It  is  interesting  to 
note  that  the  parameters  at  e/3p  =  0.6  match  those  for  the  ST-IB  reference  case, 
demonstrating  a  linkage  between  ST-IB  and  ST-2B  as  ejSp  varies  from  0.6  to  1.0. 

The  linkage  between  the  two  reference  caises  indicates  an  operating  path  may 
exist  between  high  beta  operation  in  the  1st  stability  regime  and  operation  in  the 
theoretical  2nd  stability  regime  with  beta  values  approaching  or  exceeding  unity. 
This  is  important  when  considering  what  the  access  path  might  be  from  ignition 
to  steady-state  operation  in  the  2nd  stability  regime.  While  reactor  performeince 
increzises  significantly  as  e/3p  increases  from  0.5  to  1.0,  higher  values  produce  only 
slightly  better  performance  overall,  as  seen  in  the  COE  curve. 

3.4  Variation  of  Density  and  Temperature  Profile 

The  nominal  values  adopted  for  the  density  profile  factor  (a„  =  0.5)  and 
the  temperature  profile  factor  (ay  =  1.0)  represent  reasonable  estimates  for  initial 
study  (see  A.5).  These  would  have  to  be  refined  through  detailed  MED  analysis 
to  ensure  that  the  plasma  pressure  profile,  defined  by  the  density  and  temperature 
profiles,  is  consistent  with  the  current  density  profile,  which  was  not  included  in  the 
model  used  for  this  study.  The  radial  plasma  pressure  is  proportional  to  density 
and  temperature; 

p{p)  (xn{p)T{p)  =  {n){T){l  +a„)(l  +  ar)(l 

Assuming  the  pressure  profile  defined  in  the  reference  case  to  be  valid,  it  must  remain 
constant  as  the  density  and  temperature  profiles  vary,  so  that  MED  equilibrium  and 
stability  is  maintained. 


38 


The  eissumed  pressure  profile  can  remain  relatively  const2mt  by  making  the 
simplified  assumption  that  a„  +  ot  remains  constant.  Figures  3.7  and  3.8  show 
the  effect  on  reactor  performance  for  q„  values  over  the  range  0. 1-1.3,  where  a-p  = 

1.5  —  a„.  These  figures  show  slight  improvement  in  reactor  performance  with  more 
peaked  density  profiles  and  flatter  temperature  profiles.  This  indicates  there  is  some 
flexibility  in  the  pressure  profile  without  a  significant  impact  on  COE.  If  technology 
cam  be  developed  to  form  and  sustain  peaked  density  profiles  (see  A.5),  this  may  be 
very  useful  for  overall  profile  control  (i.e.,  selecting  the  best  combination  of  density 
and  temperature  profiles  for  optimum  MHD  equilibrium  and  stability). 

3.5  Variation  of  the  Fuel  Ratio 

The  nominal  value  for  the  fuel  ratio,  na/nxj  =  1.0,  represents  an  apparent 
optimum  for  D-^He  fusion.  However,  since  the  maximum  fusion  power  density  is 
virtually  independent  of  the  fuel  ratio  (see  A.6),  reactor  performance  sensitivity  to 
variations  may  be  oriented  towards  other  important  design  parameters.  Figures  3.9 
and  3.10  show  the  effect  on  reactor  performance  for  nsfnu  values  over  the  reinge 
0. 5-2.0.  The  neutron  wall  load  curve  clearly  indicates  the  lower  irradiation  levels 
associated  with  ®He-rich  fuel  ratios.  It  is  worth  noting  that  the  ®He/D  injection  rate 
ratio  is  almost  identical  to  the  effective  fuel  ratio  in  all  fuel  cycles  considered  here 
(i.e.,  these  plasmas  are  loss  dominated  rather  than  fusion  dominated  [4]). 

These  figures  indicate  reactor  performance  is  minimized  for  fuel  ratios  around 
0.75,  with  improved  performance  for  both  higher  and  lower  fuel  ratios,  but  do  not 
address  plasma  ignition.  Increase  of  the  fuel  ratio  much  beyond  2.0  reduces  the 
neutron  wall  load  somewhat,  but  leads  out  of  the  ignition  domain  [4].  It  is  worth 
noting  here  that  the  shielding  requirement  does  not  reduce  Hnearly,  but  roughly 
logarithmically,  with  the  neutron  wall  loading.  Decrease  of  the  fuel  ratio  beyond 
0.5  simply  increases  the  neutron  wall  load  to  unacceptably  high  levels,  regardless  of 
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ignition  requiiements. 

3.6  Variation  of  Confinement  Time  Ratio,  tJte 

The  nominal  value  for  the  ratio  of  peirticle  confinement  time  to  energy  confine¬ 
ment  time,  Tp/r^  3.0,  represents  a  reasonable  estimate  for  initial  study  (see  A.6). 
(ARIES-III,  by  comparison,  assumes  Tp/te  ^  2.0).  Actual  particle  transport  phe¬ 
nomena  in  any  given  reactor  configuration  could  be  different,  and  means  to  actively 
manipulate  conductive  power  losses  are  not  well  understood.  Figure  3.11  shows  the 
effect  on  reactor  performance  for  Tp/r^  values  over  the  range  1.5-5. 0.  These  figures 
indicate  only  minor  effects  on  performance  over  the  expected  range  2  <  Tp/r^  <  4. 
This  confirms  the  expectation  that  ST  reactors  operating  at  high  betr.  Eire  more 
robust  to  uncertainties  in  particle  confinement  and  Eish  buildup  effects,  as  well  as 
synchrotron  radiation  effects  [19].  This  is  related  to  the  relatively  lower  radiation 
power  losses  at  higher  beta  levels,  which  reduces  the  impact  of  ash  buildup. 

3.7  Variation  of  Confinement  H-factor 

The  1990  ITER  H-mode  pow^  scaling  used  for  this  study  determines  the 
confinement  time  enhancement  factor  required  for  global  power  balance  (see  A.9). 
The  required  H-factor  for  ST-2B  was  3.633,  which  represents  an  optimistic  value. 
Actual  confinement  may  be  better,  but  worse  confinement  should  be  anticipated. 
Figures  3.12  and  3.13  show  the  effect  on  reactor  performance  for  decreasing  H-factor 
values  over  the  range  3.63-3.0.  These  figures  show  a  sharp  dechne  in  performance 
as  confinement  worsens  and  erratic  behavior  for  values  below  3.3.  Injection  power 
requirements  and  TF  coil  resistive  losses  both  increase  for  operation  at  lower  H- 
factors,  which  drives  up  the  recirculating  power  fraction  and  COE. 
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3.8  Variation  of  Bootstrap  Current  Fraction 

The  assumed  bootstrap  current  fraction,  Ibtl^p  =  0.90,  represents  an  opti¬ 
mistic  meudmum  value.  At  least  5-10%  of  the  plasma  current  should  be  extemeilly 
driven  to  provide  the  seed  current  for  the  bootstrap  effect  and  to  provide  current 
density  profile  control.  The  bootstrap  current  fraction  for  the  ST-2B  configuration 
may  be  less  than  90%,  based  upon  its  actual  pressure  profile  and  paxtide  diffusion 
effects.  Figures  3.14  eind  3.15  show  the  effect  on  reactor  performance  for  decreasing 
values  of  the  bootstrap  current  fraction  over  the  range  0.90-0.60.  These  figures  show 
a  Hnear  decline  in  performance  proportional  to  the  bootstrap  current  fraction.  As 
more  power  is  required  for  current  drive,  the  recirculating  power  fraction  and  COE 
increase.  The  plasma  temperature  also  increases,  and  possible  limitations  on  tem¬ 
perature  (or  temperature  gradient)  may  impose  a  minimum  acceptable  bootstrap 
current  fraction  necessary  for  safe  operation.  If  a  peaked  density  profile  and  a  rel¬ 
atively  flat  temperature  profile  can  be  formed,  as  previously  discussed,  then  high 
bootstrap  current  fractions  (80-90%)  at  acceptable  operating  temperatures  might 
be  possible. 

3.9  Variation  of  Aspect  Ratio 

The  aspect  ratio  lower  limit  of  1.2  adopted  for  the  study  represents  an  arti¬ 
ficially  imposed  bound,  allowing  enough  room  for  an  adequate  centerpost  and  to 
accommodate  for  engineering  difficulties.  Theoretically,  this  Hmit  could  be  lower. 
However,  unforeseen  physics  and  engineering  constraints  coidd  impose  a  higher 
bound.  Figures  3.15  and  3.16  show  the  effect  on  reactor  performance  for  aspect 
ratios  over  the  range  1.16-1.30. 

The  most  significant  effect  indicated  by  these  figures  is  the  linear  increase  in 
magnetic  field  strength  and  TF  coil  resistive  power  proportionEd  to  aspect  ratio.  The 
COE  curve  generally  shows  a  linear  increase,  but  the  total  capital  cost  curve  shows 
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only  a  marginal  decrease  at  lower  aspect  ratio  and  a  linear  increase  with  higher 
aspect  ratio  values.  The  plasma  beta  curve  reflects  the  higher  beta  levels  at  lower 
aspect  ratio,  as  usually  observed,  but  shows  a  leveUng  off  at  aspect  ratios  below 
1.2.  The  maximum  current  density  in  the  centerpost  tends  to  limit  any  increases  in 
reactor  performance  at  these  very  low  aspect  ratios,  but  is  not  so  critical  at  higher 
aspect  ratios. 

3.10  Variation  of  Centerpost  Current  Density 

The  maximum  current  density  in  a  single  dispersion  strengthened  copper 
conductor  centerpost  is  estimated  to  be  75  MA/m*,  which  was  the  value  adopted 
for  the  study.  Whether  such  a  component  can  be  manufactured  to  this  spedflcation 
or  greater,  capable  of  withstanding  high  stresses  and  neutron  irradiation  over  its 
lifetime,  is  uncertain.  Even  if  the  technological  uncertainties  can  be  resolved,  such  a 
component  might  not  be  cost  effective  rdative  to  a  centerpost  with  a  lower  current 
density. 

Figures  3.18  and  3.19  show  the  effect  on  reactor  performance  for  centerpost 
current  density  values  over  the  range  65-85  MA/m^.  These  figures  generally  show  a 
hnear  improvement  in  performance  proportional  to  the  current  density.  If  current 
densities  less  than  75  MA/m*  are  required,  the  impact  on  COE  is  not  very  significant. 
The  increases  in  reactor  performance  for  current  densities  greater  then  75  MA/m* 
eire  minimal,  which  is  probably  due  to  the  aspect  ratio  lower  bound  of  1.2.  High 
centerpost  current  densities  and  very  low  aspect  ratios  are  inextricably  linked,  as 
also  observed  in  the  previous  section.  Further  studies  should  consider  these  two 
parameters  together  for  more  meaningful  results. 
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3.11  Variation  of  Thermal  Conversion  Efficiency 

The  40%  thermaJ  conversion  efficiency  assumed  for  the  study  represents  a 
reasonable  value,  based  upon  industry  experience.  The  thermsd  conversion  efficiency 
is  an  average  based  upon  the  heat  energy  removal  from  the  shield  eind  divertor 
plates.  The  efficiency  could  be  higher  in  the  shield  (ARIES-III  zissumes  44%),  but 
high  efficiency  in  the  divertor  plates  might  not  be  attainable.  Figure  3.20  shows 
the  effect  on  reactor  performance  for  average  thermal  conversion  efficiency  values 
over  the  range  30-45%.  These  figures  show  a  linear  improvement  in  performance 
proportional  to  efficiency,  which  is  to  be  expected.  Actually,  this  vziriation  shows  a 
larger  impact  on  COE  than  many  of  the  "physics  uncertainty”  variations  discussed 
previously. 

3.12  Variation  of  Direct  Energy  Conversion  Efficiency 

A  preliminary  estimate  of  the  efficiency  of  a  direct  electrodynamic  convertor 
is  50%,  which  was  the  value  adopted  for  the  study.  Details  of  the  particle  drifts 
and  distribution  functions  within  the  DEC  need  to  be  resolved,  but  the  basic  design 
appears  feasible.  Figure  3.21  shows  the  effect  on  reactor  performance  for  DEC  effi¬ 
ciency  values  over  the  range  40-60%.  These  figures  also  show  a  linear  improvement 
in  performance  proportional  to  efficiency.  However,  comparing  the  COE  curves 
in  Figures  3.20  and  3.21,  the  sensitivity  to  DEC  efficiency  is  not  as  great  as  that 
to  thermal  efficiency.  Thus,  a  conceptual  direct  electrodynamic  convertor  with  an 
efficiency  of  only  40%  would  not  degrade  reactor  performance  significantly. 

3.13  Summary 

These  parametric  investigations  have  shown  the  sensitivity  of  reactor  perfor¬ 
mance  to  several  key  design  parameters.  The  ST  fusion  reactor  appears  to  be  very 
sensitive  to  the  edge  safety  factor  hmit,  the  hmit,  the  confinement  H-factor,  and 
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the  bootstrap  current  fraction.  It  appears  to  be  more  robust  with  regard  to  the 
other  parameters  studied,  permitting  considerable  flexibility  without  significantly 
affecting  reactor  performance. 
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Parameter 

ST-IA 

ST-2A 

ARIES-m 

Aspect  ratio,  A  =  RqI a 

1.20 

1.20 

3.00 

Major  radius,  Ro  (m) 

2.676 

2.344 

7.500 

Minor  radius,  a  (m) 

2.230 

1.953 

2.500 

Plasma  vertical  elongation,  kx  =  b/a 

2.665 

2.665 

1.842 

Plasma  triangularity,  6x 

0.766 

0.766 

0.814 

Edge  safety  factor,  q 

7.681 

7.681 

6.850 

Troyon  coefficient,  gj 

2.84 

4.74 

15.10 

Plasma  toroidal  beta,  0 

0.593 

0.988 

0.239 

0.60 

1.00 

1.803 

Electron  n-temperature,  [Tg]  (keV) 

66.12 

71.10 

53.29 

Average  electron  density,  (ne)(10“/m®) 

4.331 

4.658 

3.272 

Average  ion  density,  (nj)(10“/m*) 

3.193 

3.463 

2.075 

Fuel  ratio,  rug^lriD 

1.00 

1.00 

1.07 

Density  profile  factor,  a„ 

0.50 

0.50 

0.140 

Temperature  profile  factor,  ax 

1.00 

1.00 

0.792 

Plasma  effective  charge, 

2.217 

2.232 

1.998 

Lawson  parameter,  nr£;(10*^s/m®) 

1.88 

1.81 

2.45 

Particle  confinement  time,  Tp  (s) 

13.03 

11.67 

22.55 

Ratio  Tp/r^; 

3.0 

3.0 

2.0 

ITER90  H-factor,  fi 

2.099 

3.014 

2.086 

On-axis  toroidal  magnetic  field.  Bo  (T) 

2.624 

2.200 

7.575 

Peak  TF-coil  magnetic  field.  Be  (T) 

18.19 

15.59 

13.971 

Plasma  toroidal  current.  Ip  (MA) 

122.0 

89.54 

29.88 

Bootstrap  current  fraction,  Ips/Ip 

0.90 

0.90 

0.911 

NBI  beam  energy,  Ebeam  (MeV) 

5.29 

4.93 

18.40 

Current  drive  efficiency,  7(10*°A/W  —  m^) 

2.336 

2.458 

1.066 

Absorbed  current-drive  power,  Pcd  (MW) 

68.10 

40.97 

116.93 

Total  fusion  power,  Pp  (MW) 

3261 

2759 

2682 

Plasma  Q- value,  Qp  =  Pp/PcD 

47.88 

67.34 

15.59 

Radiation  power  loss,  P^ad  (MW) 

1835 

1490 

1967 

Neutron  power,  Pn  (MW) 

210 

177 

109 

Neutron  wall  loading  {MW/m^) 

0.441 

0.483 

0.079 

Gross  electric  power  (MWe) 

1844 

1571 

1319 

Net  electric  power  output  (MWe) 

1000 

1000 

1000 

Recirculating  power  fraction 

0.460 

0.360 

0.242 

Thermal  conversion  efficiency,  Tjth 

0.40 

0.40 

0.44 

Magnet  costs  ($M) 

668 

368 

253 

Current  drive  costs  ($M) 

380 

244 

478 

First- wall/ shield  costs  ($M) 

138 

108 

181 

Total  cost  ($B) 

7.20 

5.30 

3.63 

Cost  of  electricity,  COE  (mill/kWh) 

116.5 

88.5 

74.3 

Table  3.1:  Key  parameters  of  the  1st  Stability  ST,  2nd  Stability  ST  and  ARIES-III 
2nd  Stability  Tckamak 


Plasma  Currenl,  Ip  (MA)  Major  Radius,  Ro  (m) 
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Edge  Sa&tj  Factor, 


Edge  Safety  Factor, 


Edge  Safety  Factor,  q^ 


Edge  Safety  Factor,  q^ 


Figure  3.1:  ST-IA  Sensitivity  to  Variation  of  Edge  Safety  Factor.  Ro,Bq,Ip,  and 
COE  versus 


PiBima  Current,  7p  (MA)  Major  lUdiui,  Ro  (m) 
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Edge  Safety  Factor,  g*  Edge  Safety  Factor,  g* 


Figure  3.3:  ST-2A  Sensitivity  to  Variation  of  Edge  Safety  Factor:  Ro,Bo,Ip,  and 
COE  versus 


Pint  (MW) 


Edge  Sainy  Factor, 


Edge  Sa&ty  Factor, 


Figure  3.4:  ST-2A  Sensitivity  to  Variation  of  Edge  Safety  Factor:  Pj„j,  TF  Coil 
Resistive  Power,  Recirculating  Power  Fraction,  and  Total  Cost  versus 
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Parameter 

ST-IB 

ST-2B 

Aspect  ratio,  A  =  Ro/d 

1.20 

1.20 

Major  radius,  Ro  (m) 

2.154 

1.912 

Minor  radius,  a  (m) 

1.795 

1.594 

Plasma  vertical  elongation,  kx  =  b/a 

2.665 

2.665 

Plasma  triangularity,  Sx 

0.766 

0.766 

Edge  safety  factor,  q 

5.0 

5.0 

Troyon  coefficient,  gy 

2.06 

3.43 

Plasma  toroidal  beta,  /3 

0.660 

1.100 

0.60 

1.00 

Electron  n-temperature,  [T,,]  (keV) 

71.85 

75.52 

Average  electron  density,  (ne)(10“/Tn®) 

5.086 

5.640 

Average  ion  density,  (nj)(10“/m®) 

3.776 

4.223 

Fuel  ratio,  nagginjj 

1.00 

1.00 

Density  profile  factor,  a„ 

0.50 

0.50 

Temperature  profile  factor,  ay 

1.00 

1.00 

Plasma  effective  charge,  Z^ff 

2.228 

2.244 

Lawson  parameter,  nr^(10*^s/m^) 

1.77 

1.76 

Particle  confinement  time,  Tp  (s) 

10.46 

9.37 

Ratio  Tj,/te 

3.0 

3.0 

ITER90  H-factor,  fi 

2.582 

3.633 

On-axis  toroidal  magnetic  field,  Bo  (T) 

1.955 

1.649 

Peak  TF-coil  magnetic  field,  Bg  (T) 

14.09 

12.19 

Plasma  toroidal  current.  Ip  (MA) 

112.4 

84.18 

Bootstrap  current  fraction,  last  Ip 

0.90 

0.90 

NBI  beam  energy,  (MeV) 

5.026 

4.995 

Current  drive  efficiency,  7(10*° A/PP  —  m*) 

2.489 

2.587 

Absorbed  current-drive  power,  Pcd  (MW) 

52.85 

36.09 

Fusion  power,  Pp  (MW) 

2585 

2348 

Pleisma  Q-value,  Qp  =  Pp/PcD 

48.92 

65.06 

Radiation  power  loss.  Prod  (MW) 

1384 

1231 

Neutron  power,  (MW) 

165 

150 

Neutron  wall  loading  (MW/m^) 

0.535 

0.616 

Gross  electric  power  (MWe) 

1488 

1353 

Net  electric  power  output  (MWe) 

1000 

1000 

Recirculating  power  fraction 

0.330 

0.260 

Thermal  conversion  efficiency,  r/th 

0.40 

0.40 

Magnet  costs  ($M) 

398 

231 

Current  drive  costs  ($M) 

304 

220 

First-wall/shield  costs  ($M) 

92 

74 

Total  cost  ($B) 

5.33 

4.29 

Cost  of  electricity,  COE  (mill/kWh) 

88.9 

73.7 

Table  3.2:  Key  parameters  of  the  1st  Stability  ST  and  2nd  Stability  ST  reference 


cases 


Platma  Bela, 


Plasma  Currcnl,  h  (MA)  Major  nmliiis,  Rq  (ni) 


necirc.  Power  FVac.  (xlO~')  Deniity,  (n,)  (10*®m 
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Density  Profile  Fsctot, 


Density  Profile  Factor,  a„ 


Density  Profile  Factor, 


Density  Profile  Factor,  On 


Figure  3.8:  ST-2B  Sensitivity  to  Variation  of  Density  and  Temperature  Profile:  {rig), 
(Te),  Recirculating  Power  Fraction,  and  Total  Cost  versus  a„  (where  a^-  —  1.5  —  a„) 


PIa.m.  Current,  h  (MA)  Major  lladiu..  (m) 


Neutron  Well  Load  (MW/m*)  Deniity,  (n,)  (10**m 
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Fad  Ratio,  nz/ng 


1.2 
1.0 
0.8 
0.6 
0.4 
0.2 

0.5  1.0  1.5  2.0 

Fad  Ratio,  nj/no 


Fad  Ratio,  nilug 


Fad  Ratio,  n^lng 


Figure  3.10:  ST-2B  Sensitivity  to  Variation  of  Effective  Fuel  Ratio:  (ne),Pjp,  Neu¬ 
tron  WzJl  Load,  and  Total  Cost  versus  na/n^j 


Current 


Coniineinent  H-£actor, 


Coafineraeat  H-factor,  /l 


Figure  3.12;  •ST-2B  Sensitivity  to  Variation  of  Confinement  H-Factor:  Ro,Bo,Ip, 
and  COE  versus  fi 
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Figure  3.13;  ST-2B  Sensitivity  to  Variation  of  Confinement  H-Factor;  Fi„j, 
Recirculating  Power  Fraction,  and  Total  Cost  versus  fi 
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Figure  3.14:  ST-2B  Sensitivity  to  Variation  of  Bootstrap  Current  Fraction: 
Rq,Bo,Ip,  and  COE  versus  I},,/ Ip 


Recirc.  Power  FVac.  (xlO'*)  Injection  Power, (MW) 


Plasma  Current,  /p  (MA)  Major  Radius,  iio  (m) 


Aspect  Ratio,  Ro/a 


Aspect  Ratio,  Ro/a 


Figure  3.16:  ST-2B  Sensitivity  to 
versus  Ro/a 
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Figure  3.17:  ST-2B  Sensitivity  to  Variation  of  Aspect  Ratio:  /3,  TF  Coil  R.esistive 
Power,  Recirculating  Power  Fraction,  and  Total  Cost  versus  Ro/a 
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Figure  3.18:  ST-2B  Sensitivity  to  Variation  of  Centerpost  Current  Density: 
i2o,  Bo, /p,  and  COE  versus 


Recirc.  Power  Frac.  (xlO“*)  Plaima  Beta, 
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Figure  3.19:  ST-2B  Sensitivity  to  Veuiation  of  Centerpost  Current  Density:  /?,  TF 
Coil  Resistive  Power,  Recirculating  Power  Fraction,  and  Total  Cost  versus  J_ 
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Figure  3.20:  ST-2B  Sensitivity  to  Variation  of  Thermal  Conversion  Efficiency; 
Ro,Bo,Ip,  and  COE  versus  7/t^ 
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Figure  3.21:  ST-2B  Sensitivity  to  Variation  of  Direct  Energy  Conversion  Effiaency 
Ro,Bo,Ip,  and  COE  versus 


CHAPTER  4 

Discussion  and  Conclusions 


4.1  Introduction 

Based  on  the  ST-2B  reference  case  used  for  the  study,  compact  (iZo  <  2m)  D- 
^He  fusion  at  low  field  (S©  <  2T)  and  modest  cost  (COE  ~  74  mitl/kWh,  comparable 
to  ARIES-III)  appears  viable.  The  assumptions  adopted  for  the  study  represent  rea¬ 
sonable  extrapolations  of  physics  and  technology,  which  might  be  attainable  by  the 
time  a  D-^He  ST  demonstration  plant  could  become  operational  {~  20507).  The 
parametric  investigations  discussed  in  Chapter  3  describe  the  reactor  performance 
sensitivities  to  several  key  parameters.  The  supporting  database  for  those  param¬ 
eters  can  be  refined  as  future  experiments  resolve  physics  issues  and  technological 
advances  occur. 

Assuming  the  ST-2B  reference  case  represents  the  anticipated  parameters  of  a 
D-®He  ST  fusion  reactor,  it  is  useful  to  aissess  its  economic,  safety  and  environmental 
characteristics.  Among  the  specific  reference  cases  eicamiued  by  the  Senior  Com¬ 
mittee  on  Environmental,  Safety,  and  Economic  Aspects  of  Magnetic  Fusion  Energy 
(ESECOM)  [1],  the  D-^He  tokamak  was  one  of  the  most  potentially  attractive,  but 
technological  uncertainties  and  lack  of  design  detail  made  its  feasibility  uncertain. 
The  ESECOM  anzilysis  emphasized  future  development  of  reactor  designs  combining 
high  levels  of  safety  assurance,  high  mass  power  density,  direct  energy  conversion, 
and  design  simphdty  for  rehabUity  and  ease  of  maintenance  necessary  to  realize 
the  possibility  of  cost-competitive  electricity.  The  D-®He  ST  meets  these  challenges 
head-on.  With  such  well-defined  criteria  available,  the  potential  advantages  of  the 
ST-2B  reference  case  may  be  assessed  realistically  and  quantitatively. 
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4.2  Economic  Characteristics 

The  COE  of  the  ST-2B  does  not  appear  to  be  economically  competitive  with 
advanced  fission  breeder  reactors,  but  does  appear  to  be  economically  competitive 
with  large-scale  solar  electric  plants,  as  discussed  in  Chapter  3  (see  3.2.1).  There  are 
many  uncertainties  involved  in  the  COE  estimate  for  ST-2B.  The  uncertainties  with 
the  largest  effects  on  ESECOM’s  estimates  of  COE  were  associated  with  achiev¬ 
ing  the  assumed  plasma  beta  and  current-drive  efficiency  and  with  attaining  plant 
rehabihty  and  maintainabihty  chciracteristics  compatible  with  reasonable  capacity 
factors.  These  are  crucial  issues  for  the  D-®He  ST,  too. 

The  design  characteristics  [1]  offering  the  most  important  potential  benefits 
for  fusion  COE  are:  (1)  compactness,  which  reduces  the  capital  cost  of  the  fusion 
power  core  and  the  sensitivity  of  COE  to  plasma  performance,  and  may  ease  mainte¬ 
nance;  (2)  high  level  of  public  safety  assurance  based  on  low  radioactive  inventories 
and  passive  mechanisms  for  preventing  releases,  which  should  reduce  costs  for  ac¬ 
tive  safety  systems  and  nucleEir-grade  components  as  weU  as  facilitating  siting  and 
hcensing;  and,  (3)  advanced  energy  conversion  systems,  which  should  be  able  to 
reduce  balance  of  plant  (BOP)  costs  and  may  increase  capacity  factors. 

The  ST-2B  reference  case  combines  all  three  features  mentioned  above  in  one 
design  (ARIES-III  is  not  compact  and  does  not  use  advanced  energy  conversion).  If 
design  simplification  C2in  reduce  the  number  of  plant  components  and  design  rugged¬ 
ness  can  enhance  critical  design  margins  and  extend  plant  hfe  (>  60  years),  then 
additional  reductions  in  the  COE  estimate  may  be  possible  for  ST-2B.  These  favor¬ 
able  features  might  be  able  to  offset  some  physics  and  technology  shortcomings,  and 
permit  a  D-®He  ST  fusion  reactor  to  remain  economically  competitive.  Additionally, 
it  might  be  possible  to  design  smaller  modrilar  reactors  (100-500  MWe),  where  the 
potenticil  for  mass  production  could  lower  the  capital  cost  per  kilowatt  capacity.  A 
low  initial  capital  cost  would  then  permit  the  spread  of  nuclear  energy  to  developing 
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countries  and  smaller  communities  of  developed  countries. 

4.3  Safety  and  Environmental  Characteristics 

The  inherent  advantages  of  the  D-®He  fuel  cycle  (i.e.,  reduced  neutron  acti¬ 
vation  and  tritium  handhng  problems  relative  to  the  D-T  fuel  cycle)  signiiic£intly 
reduces  the  safety  eind  environmental  risks  normally  associated  with  D-T  fusion  and 
advanced  fission  reactors.  ESECOM’s  analysis  of  safety  and  environmental  char¬ 
acteristics  assessed  the  following  areas:  (1)  possibihties  and  consequences  of  major 
releases  of  radioactivity  from  reactor  accidents;  (2)  magnitude  of  the  radioactive 
waste  burden;  (3)  occupational  and  public  exposures  to  radiation  in  routine  oper¬ 
ation;  and,  (4)  unwanted  links  to  nuclear  weaponry.  The  next  few  sections  assess 
these  characteristics  qualitatively  and  quantitatively  for  the  ST-2B  reference  caise. 

4.3.1  Risk  of  Major  Releases  of  Radioactivity 

The  probability  and  consequences  of  major  radioactivity  releases  during  ac¬ 
cidents  for  any  nuclear  power  plant  must  be  covered  by  design  and  demonstrated 
during  the  licensing  procedure.  The  largest  quantities  of  radioactivity  are  generally 
in  the  form  of  neutron  activation  products,  most  of  which  are  embedded  in  sohd  re¬ 
actor  materials.  Since  the  neutron  wall  load  for  ST-2B  is  about  0.6  MW/m^,  which 
is  about  25%  that  of  a  D-T  tokamak,  activation  products  in  the  reactor  structure 
would  accumulate  more  slowly.  The  steady-state  inventory  of  radioactive  tritium  in 
the  ST-2B  plasma  is  less  than  1  g  and  the  total  inventory  in  the  tritium  handhng 
system  is  probably  much  less  than  1  kg,  but  this  is  notably  less  than  the  anticipated 
tritium  inventories  for  D-T  reactors  (up  to  5  kg).  Such  a  low  tritium  inventory  in 
ST-2B  represents  a  radioactive  hazard  an  order  of  magnitude  less  than  a  D-T  reac¬ 
tor  and  many  orders  of  magnitude  less  than  the  fuel  inventory  hazard  in  a  fission 


reactor. 
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Prominent  features  of  fusion  in  general  with  regard  to  safety  are  the  expected 
self-termination  of  the  fusion  process  under  accidental  conditions  and  the  fact  that 
a  recriticality  problem  does  not  exist.  Nevertheless,  the  critical  safety  issues  of 
concern  regarding  fusion  reactors  axe  [31]:  (1)  loss  of  cooling  accidents;  (2)  loss 
of  vacuum  accidents;  (3)  tritium  system  failure;  and,  (4)  magnet  system  failure. 
Without  elaborating  on  the  details  of  these  issues,  if  an  accident  in  ST-2B  occurs,  the 
total  radioactivity  release  to  the  environment  in  the  form  of  tritium  and  activation 
products  can  be  kept  below  regulatory  limits  [31]. 

The  ST-2B  reference  case  was  assumed  to  have  a  level  of  safety  assurance 
(LSA)  of  2.  In  a  level  2  reactor,  safety  is  assured  by  passive  mechanisms  of  release 
limitation  as  long  as  severe  reconfiguration  of  large-scale  geometry  is  avoided,  and 
escalation  to  fatality-producing  reconfigurations  from  less  severe  initiating  events 
can  plausibly  be  precluded  by  passive  design  features  [1],  The  quantitative  threshold 
for  avoidance  of  early  fatalities  off-site  is  a  critical  whole-body  dose  of  200  rem  at 
the  plant  boundary,  assumed  to  be  1  km  from  the  point  of  the  release.  Level  2 
is  probably  a  conservative  assumption,  considering  the  low  level  of  radioactivity  in 
ST-2B.  It  may  be  possible  to  warrant  the  classification  of  ST-2B  as  Level  1,  where 
safety  is  assured  by  passive  mechanisms  of  release  limitation  no  matter  what  the 
accident  sequence.  For  either  classification,  however,  the  passive  safety  features  and 
the  reduced  probability  of  a  severe  accident  (both  in  terms  of  public  safety  and 
investment  protection)  should  significantly  ease  the  siting  and  bcensing  process  for 
a  D-®He  ST  power  plant. 

4.3.2  Radioactive  Waste  Burden 

The  only  radioactive  wastes  associated  with  ST-2B  are  those  components 
which  will  contain  neutron  activation  products,  primarily  the  center  conducting 
post  and  the  first  wall/shield.  The  fully  exposed  centerpost  will  receive  neutron 
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damage  and  require  replacement  at  regular  intervals,  perhaps  every  few  years  [18]. 
The  activated  centerpost,  however,  contains  no  sensitive  or  toxic  fissile  material, 
requires  httle  processing,  and  is  self-contained  *i  a  moaast  volume  (~  10.3m®).  De¬ 
spite  the  low  neutron  wall  load,  the  first  wall  and  shield  will  probably  accumulate 
a  high  level  of  radioactivity  over  their  effective  lifetin>es,  which  will  be  longer  for 
ST-2B  than  for  D-T  reactors.  Two  separate  studies  [24,  32]  indicate  that  if  low 
activation  components  are  used  for  the  first  wall  and  shield,  the  criteria  for  Class 
A  disposal  according  to  10CFR61  would  be  satisfied.  Class  A  refers  to  minimal  re¬ 
striction  surface  storage.  Components  made  of  conventional  materials  such  as  Type 
316  stainless  steel  and  Fel422  steel  would  require  Class  C  (restricted)  surface  stor¬ 
age.  In  either  case,  these  wastes  would  qualify  for  shallow  burial  under  the  current 
regulatory  philosophy.  There  would  be  no  high-level  waste,  such  as  the  spent  fuel 
elements  from  fission  reactors,  that  requires  deep  geological  disposal.  The  overall 
reduced  level  and  volume  of  radioactive  waste  for  ST-2B  would  permit  on-site  dis¬ 
posal  until  deactivation  to  acceptable  levels,  presenting  minimal  risk  to  the  pubUc 
and  environment  at  minimal  cost. 

4.3.3  Exposures  to  Radiation 

The  normal  operation  of  nuclear  power  plants  entail  small  emissions  of  ra¬ 
dioactivity  to  the  environment  and  some  exposure  of  the  pubhc  to  radiation,  as 
well  as  routine  exposure  of  workers  to  radiation  at  plant  sites.  For  ST-2B,  radioac¬ 
tivity  release  to  the  environment  in  the  form  of  tritium,  ’^C,  and  other  activation 
products  would  be  neghgible  during  normal  operation,  well  below  regulator}"  hmits. 
Radiological  exposure  of  workers  normally  reaches  its  highest  levels  during  reactor 
maintenance  and  repair.  With  higher  capacity  factors  and  reduced  activation  of 
components,  radiological  exposure  of  workers  should  be  significantly  lower  than  for 
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advanced  fission  and  D-T  fusion  reactors.  If  low  activation  SiC  composites  or  Al- 
Mg-Si  alloys  with  impurities  <  1  ppm  could  be  used  as  structural  materials,  then 
limited  contact  maintenance  outside  the  shield  may  be  permissible  [33].  Unfortu¬ 
nately,  the  materials  community  cannot  identify  with  confidence  any  ceramic  or 
alloy  that  meets  all  the  diverse  requirements  as  a  structural  material  for  a  fusion 
reactor  (i.e.,  fabricability,  mechanical  strength,  chemical  compatibibty,  cost,  etc.) 
[34]. 


4.3.4  Nuclear  Weapons  Link 

In  general,  fusion  energy  systems  are  less  likely  than  fission  energy  systems  to 
contribute  to  the  acquisition  of  nuclear  weapons  capabilities  by  subnational  groups 
because  is  not  a  limiting  constraint  on  thermonuclear  weapon  construction 

the  way  fissile  materials  are  a  constreiint  on  fission  weapon  construction  [1].  Since 
the  tritium  inventory  in  ST-2B  is  very  small,  even  this  link  to  nuclear  weapons  is 
tenuous.  In  contrast,  the  relative  weapons  potential  of  an  advanced  fission  breeder 
reactor  is  very  high  and  represents  a  risk  to  international  security.  It  may  be  ad¬ 
vantageous,  however,  to  allow  tritium  to  decay  to  ^He  before  it  is  reinjected  into 
the  plasma  (as  assumed  for  this  study)  at  the  expense  of  tritium  storage  require¬ 
ments  which  safeguard  against  weapons  proliferation.  In  general,  though,  a  D-^He 
ST  fusion  reactor  presents  minimal  risk  and  may  be  freely  promoted  throughout  the 
world  as  a  relatively  safe  power  source  with  an  unhmited  fuel  supply. 

4.4  Conclusion 

The  anticipated  DAHe  spheric^ll  torus  fusion  reactor  demonstrates  attractive 
economic,  safety  and  environmental  characteristics  which  should  justify  its  continued 
research  and  development  (Ri:D).  This  conclusion  is  based  on  assumptions  about 
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plasma  performance  and  engineering  characteristics  that  are  optimistic  but  defensi¬ 
ble  extrapolations  from  current  experience  and  apphcation  of  a  computer  model  to 
quantitatively  assess  the  performance  of  a  reference  case  configuration.  Overall,  the 
D-^He  ST  reactor  appears  to  be  competitive  with  advanced  fission  breeder  reactors 
and  large-scale  solar  electric  plants  as  a  future  energy  source  for  the  long  term. 

The  physics  and  engineering  extrapolations  adopted  for  this  system  study 
highhght  a  number  of  outstanding  issues  which  need  to  be  resolved  through  con¬ 
tinued  R&D  of  ST  reactors.  These  include:  (1)  tokamak-like  confinement  at  very 
low  aspect  ratio;  (2)  stable  of  order  unity;  (3)  plasma  stabiHty  at  low  edge  q; 
(4)  bootstrap  current  fraction  of  order  unity;  (5)  power  and  particle  control  using 
expanded  poloidal  divertors  vith  gaseous  targets;  (6)  particle  drifts  and  distribution 
functions  within  a  direct  electrodynamic  convertor;  (7)  high  energy  neutral  beam 
current  drive  at  high  plasma  densities;  (8)  high-speed  fuel  pellet  injection  at  high 
plasma  densities  and  temperatures;  (9)  low- activation  structural  materials;  and,  (10) 
®He  mining  on  the  moon.  These  issues  represent  significant  challenges  which  must 
be  overcome  before  a  D-^He  ST  reactor  can  be  designed  with  confidence. 

Two  other  aspects  of  D-^He  ST  reactor  operation  which  will  have  to  be  ad¬ 
dressed  are  plasma  ignition  and  response  to  transients.  This  study  only  assumed 
steady-state  operation  at  high  plasma  temperature.  If  reaching  the  steady-state 
condition  requires  a  D-T  ignition  plasma,  with  its  adverse  neutron  activation  and 
tritium  handling  problems,  then  the  inherent  advantage  of  the  D-®He  fuel  cycle  wiU 
be  compromised.  If  the  D-^He  plasma  can  be  heated  directly  to  the  ignition  regime 
with  auxiliary  power  only,  perhaps  with  a  deuterium-rich  ignition  plasma  [28],  then 
D-T  assisted  start-up  could  be  avoided.  Operational  transients,  due  either  to  minor 
plasma  disruptions  or  changing  load  demands  from  the  electric  utility  grid,  require 
flexible  response  from  the  reactor  configuration.  R<kD  efforts  must  resolve  what  the 
design  margins  should  be  relative  to  the  optimum  steady-state  configuration  and 
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their  impacts  on  COE.  Since  the  reactor  performance  is  highly  non-Unear,  R&D 
efforts  must  also  resolve  how  operators  wiU  be  able  to  monitor  and  control  reactor 
performance  during  transient  scenarios.  These  are  complex  issues  that  may  only  be 
resolved  through  operational  experience  from  large  steady-state  experiments. 


APPENDIX  A 

Spherical  Torus  Physics  Model 


A.l  Introduction 

The  Spherical  TOrus  Reactor  Analysis  Code  (STORAC)  uses  a  profile- averaged 
global  model,  which  closely  follows  the  method  for  zero-dimensional  International 
Thermonuclear  Reactor  (ITER)  scoping  studies  [35],  with  the  exception  of  the  par¬ 
ticle  transport  effect  on  energy  confinement.  The  model  described  here  is  the  one 
used  for  this  study,  which  incorporates  the  chcinges  identified  in  Chapter  2.  Meter- 
kilogram- second  units  are  used  throughout  with  the  exception  of  kilo-electron-volts 
for  plasma  temperatures  and  megawatts  per  cubic  meter  for  power  densities,  unless 
otherwise  noted. 

A.2  Plasma  Shape  and  Geometry 

Free-boundary  magnetohydrodynamic  (MHD)  equilibrium  calculations  show 
that  as  the  aspect  ratio  is  reduced  to  around  1.5  for  spherical  torus  plasmas,  plasma 
elongation  of  k  2  occurs  naturally  without  a  significant  shaping  field  [16].  En¬ 
hanced  elongation  (relative  to  that  of  tolmmaks  with  conventional  aspect  ratios 
around  3)  is  a  natural  feature  of  the  spherical  torus.  MHD  calculations  also  show 
that  shaping  the  plasma  by  giving  it  a  triangularity  6  produces  a  stabilizing  effect. 
The  STORAC  scalings  for  the  plasma  sepuatrix  elongation  and  triangularity  given 
by  [19]  are: 

K  =  2.1(1.0  +  0.44e^-^) 

^  =  0.53(1.0  +  O.TTe®-®) 

where  e  is  the  inverse  aspect  ratio,  a/Ro-  This  represents  about  25%  higher  elon¬ 
gation  than  the  natural  elongation,  but  within  acceptable  hmits.  If  the  plasma  is 
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further  elongated  above  a  critical  elongation  it  becomes  unstable  to  a  global,  mainly 
vertical  instability  [11],  For  a  spherical  torus  reactor  with  aspect  ratio  around  1.2, 
e  =  0.71,  K  —  2.73  and  6  =  0.77.  The  elongation  scaling  used  here  also  fits  for  an 
ARIES-I  reactor  and  the  ITER. 

The  plasma  surface  area  Ap,  approximated  by  a  revolution  of  two  arcs  which 
intersect  the  plasma  x-points  (assuming  a  double-null  divertor  configuration)  and 
the  plasma  midplane  outer  radius  and  inner  radius,  is: 

Ap  =  47r[r<.o{(R  -I-  a  -  r^)sin~^{blr^)  4-  6}  -f  r^{{R  -  a  +  r^)sin~^{b/ra)  -  b}] 
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A. 3  Plasma  Current 


The  plasma  current  {Ip)  scaling  used  for  a  spherical  torus  is  a  function  of  the 
edge  safety  factor  (g^),  toroidal  magnetic  field  on-axis  (Bq),  and  plasma  geometry. 
In  any  toroidal  plasma  with  a  finite  aspect  ratio,  the  exact  relationship  among  q^, 
Ip,  Bq,  R,  and  a  is  complex.  For  A  <  k^/{1  ■+  S)  +  6,  which  is  normally  the  case  for 
spherical  torus  plasmas,  the  plasma  current  scahng  given  by  [19]  is 

Jp(MA)  =  ^ 

^  TT  \  E2  } 

where  the  relationship  between  the  average  poloidal  field  at  the  plasma  edge  (Rp) 
and  the  toroidal  magnetic  field  on  axis  {Bo)  is 

■^1  =  /i  -  h'lln  j  j  ^2  =  /2b  +  2h2ian~\y2)] 

f  —  <^1(1  +  ^)e  ,  _  ^2(1  —  6)€ 

^  (l+e)(cic-l)  *  (1  -  e)(c2e -f  1) 


g  =  e/c/(l  -  e^) 


1  +  (1  ~  ci)c/2 
/(I  +e)(cie  -  1) 
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1  -f  e  K 
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\/(l  -e)(c2e  f  1) 

/c^T+T  1—5 

/  1  -e  'ir~ 

2k 

d2{l  -5) 


ci=«V(l+5)  +  5  C2  =  «V(1-5)-5 

=  (ifi)  =  (r^)  +  ^ 

The  edge  safety  factor  is  an  important  parameter  which  is  related  to  the  fluid 
stability  of  the  plasma-field  configuration.  It  is  defined  as  =  aBo/ RoBp,  which 
represents  the  number  of  rotations  a  magnetic  field  hne  makes  in  the  toroidal  direc¬ 


tion  per  rotation  in  the  poloidal  rotation  (i.e.,  the  hehcal  twdsting  of  the  field  hnes 
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around  the  torus).  For  a  given  machine,  is  a  decreasing  function  of  the  plasma 
current  and  is  a  convenient  parameter  to  use  because  it  presents  a  lower  hmit  below 
which  the  MHD  stabihty  of  the  plasma  is  compromised.  Safe  operation  for  tokamaks 
is  normally  obtained  for  9,/,  >  3.0  [36].  The  default  scaling  assumed  in  STORAC  for 
the  edge  safety  factor  hmit  given  by  [19]  is 

9^  =  3.0(1  +  2.6e=®) 

For  an  aspect  ratio  of  1.2,  9^  >  7.68.  The  parameter  9^  itself  is  not  derived  in  the 
STORAC  calculations;  it  is  either  input  as  a  fixed  value  or  iterated  as  a  program 
variable  to  be  greater  than  or  equal  to  the  assumed  9^  hmit. 

It  is  also  desirable  to  have  a  simple  formula  for  the  safety  factor  which  can 
be  used  to  calculate  the  energy  confinement  time  for  various  scahngs.  The  edge 
safety  factor  is  typically  evaluated  on  the  surface  which  contains  95%  of  the  total 
flux  enclosed  by  the  separatrix.  For  a  spherical  torus,  the  assumed  MHD  q- value  at 
95%  flux  is 

9v-(95%)  1.3(1. 0  - 

For  an  aspect  ratio  around  1.2,  9,^(95%)  ~  3.41.  Plasma  confinement  tends  to 
deteriorate  with  low-q  operation  [35],  and  a  value  of  9^(95%)  —9(0)  larger  than  2  is 
assumed  to  provide  sufficient  global  shear  for  stabihty  of  kink  modes  in  ITER-hke 
reactors.  Since  the  scahngs  used  here  are  aspect  ratio  dependent,  and  may  not  be 
appropriate  at  very  low  aspect  ratio,  low-q  confinement  effects  in  a  spherical  torus 
are  unknown.  Optimistically,  lower  values  of  9^  and  9^(95%)  —  9(0)  sufficient  for 
stabihty  might  be  possible. 

A.4  Bootstrap  Current 

It  is  predicted  theoreticaUy  [37]  that  radial  diffusion  induces  a  current  in  the 
toroidal  direction  and  that  the  current  becomes  large  in  the  trapped  particle  re¬ 
gion.  In  a  plasma  with  a  finite  pressure  gradient  and  an  apphed  magnetic  field,  the 
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superposition  of  the  Lannor  orbits  of  the  plasma  particles  creates  a  diamagnetic 
’bootstrap’  current  perpendicular  to  both  the  magnetic  field  and  the  pressure  gra¬ 
dient.  For  an  axially  S3rmmetric  and  collision  free  system,  the  steady-state  current 
density  due  to  the  net  electron  flow  derived  in  [38]  is 

B,Sr 

where  £  is  the  inverse  aspect  ratio.  In  a  tokamak  plasma,  the  radial  pressure  gradient 
generates  a  substantial  current  parallel  to  the  hehcal  magnetic  field.  A  large  boot¬ 
strap  current  is  also  expected  to  be  generated  in  the  core  plasma  with  central  pellet 
deposition  and  heating,  and  in  the  edge  plasma  in  the  H-regime  [11].  A  completely 
bootstrapped  device  is  not  possible,  however,  and  a  small  ’’seed”  current  must  still 
be  driven  in  the  center  of  the  plasma,  where  the  pressure  gradient  is  neghgible.  The 
seed  current  is  then  ampHfied  by  the  bootstrap  effect  in  the  plasma. 

An  empirical  fit  for  the  bootstrap  current  fraction  assumed  for  ITER  [35]  is 

hji  = 


where 


and 


Cb.  =  1.32  -0.235[9^(95%)/g(0)]  +  0.0185[g^(95%)/g(0)]  = 

I3„  = 

B 

"  5(a) 


(a)  =  [V/i2yr^Ro)r/^ 


is  the  effective  plasma  radius. 

This  empirical  scaling  is  probably  inappropriate  for  very  low  aspect  ratios 
(relative  to  ITER).  Therefore,  the  maximum  bootstrap  current  fraction  is  optimisti¬ 
cally  assumed  to  be  90%  of  the  total  plasma  current  (ARIES-III  eissumes  91.1%). 
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It  may  be  possible  to  achieve  such  a  large  bootstrap  contribution  by  varying  the 
various  plasma  profiles  and  geometries,  as  indicated  by  the  pressure  gradient  and  Ip 
dependence  shown  in  the  previous  equations.  To  benefit  substantially  from  boot¬ 
strap  current  generation,  a  peaked  density  profile  provided  by  central  fuehng  of  the 
plasma  may  be  desirable. 


A. 5  Profile  Effects 


STORAC  uses  a  profile- averaged  global  model  to  evaluate  tokamak  perfor¬ 
mance.  A  zero- dimensional  analysis  considers  only  relationships  between  averaged 
quantities  and  ignores  any  dependence  on  poloidaJ  coordinates.  Poloidal  sections 
of  constant  density  and  temperature  are  represented  by  a  family  of  concentric  and 
similar  eUipses  given  by  [36]  as 


where  a  and  b  are  the  horizontal  and  vertical  minor  radii  of  the  plasma  poloidal 
cross  section,  respectively,  and  p  is  a  non  dimensional  radial  coordinate  equal  to  1 
at  the  plasma  boundary.  The  radial  plasma  profiles  of  density  and  temperature  are 
assumed  to  be  of  a  parabohc  form 


X{p)  =  Xo 


where  X  =  n,  T.  The  average  density,  average  temperature  and  density-weighted 
temperature  are 


{n)=lJ^ndV 


Uq 

1  +a„ 
To 

1  -i-  OCX 


^  l-fQ„ 

Jy  nd\  1  -r  a„  -h  a  J 


and 
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where  no  and  To  are  the  peaked  density  and  temperature  on  the  magnetic  axis. 
The  density  and  temperature  profiles  may  then  be  modeled  by  expressions  of  the 
following  nondimensional  form: 

n{p)  =  no(l  -  ”  =  (n)(l  +  a„)(l  -  ” 

T(p)  =  To(l  -  p^r ^  =  (T)(l  +  ar)(l  -  P^)“^ 

The  parameters  On  and  ay  represent  parabohc  profile  exponents,  where  = 
0  depicts  a  perfectly  flat  profile  and  larger  viilues  depict  increasingly  peaked  parabohc 
profiles.  The  nominal  values  a„  =  0.5  and  =  1.0,  comparable  with  ITER  assump¬ 
tions,  cire  used  as  default  values.  Actual  profiles  would  have  to  be  refined  based  on 
detailed  MHD  equihbrium  and  stabihty  analysis  to  ensure  that  the  pressure  profile 
is  consistent  with  the  current  density  profile,  which  is  not  considered  in  STORAC. 

The  relatively  fiat  density  profile  is  consistent  with  expectations  for  operation 
in  the  H-mode  regime,  where  density  profiles  tend  to  be  very  flat  and  even  inverted 
[39].  Disregarding  inverted  profiles,  0  <  <  0.75  may  be  taken  as  a  reasonable 

range.  However,  recent  experience  with  frozen  fuel  pellets  and  neutral  beam  in¬ 
jection  indicates  very  peaked  density  profiles  in  H-mode  are  possible.  A  significant 
problem  with  producing  a  fuel  pellet  conteiining  ^He  is  that  it  does  not  sohdify  unless 
it  is  pressurized  to  about  30  atmospheres.  It  might  be  possible  to  encapsulate  hquid 
^He  in  a  thin-walled  polymetric  shell  overcoated  with  frozen  deuterium,  and  inject 
such  a  peUet  at  velocities  high  enough  to  compensate  for  the  increased  ablation  rate 
associated  with  high  D-^He  plasma  temperatures  [5].  Assuming  a  peaked  density 
distribution  can  be  formed  and  heated  effectively,  peaking  factors  of  no/  (n)  =  2.5  to 
4  are  possible,  corresponding  to  1.5  <  <  3  as  a  reasonable  range.  To  sustain 

peaked  density  profiles  given  present  pellet  injector  technology,  however,  a  means  of 
creating  faster  pellets  will  have  to  be  developed. 

The  peaked  temperature  profile  appears  to  be  consistent  with  theoretical  ex¬ 
pectations.  The  temperature  profile  exponent,  qj-,  may  be  derived  in  the  following 
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way  [39].  First,  in  resistive  steady-state  with  classical  resistivity,  J  a  So, 

Jol {J)  =  1  -1-  aj  =  1  4-  la^.  This  then  implies  q^(95%)/g(0)  =  1  -f  laj.  Taking 
g(0)  =  1.0  (e.g.,  with  sawtooth  oscillations),  and  q^{95%)  ~  3.4,  we  have  Qj-  ~  1.6. 
For  the  second  stability  regime,  assuming  ^(0)  =  2.0  and  q^{9h%)  ~  3.4,  we  have 
aj  —  0.5.  This  derivation  is  not  necessaxily  true  for  beam  driven  cases,  as  assumed 
for  this  study,  so  some  corrections  must  be  allowed.  The  temperature  profile  could 
be  more  peaked  than  this,  e.g.,  with  sawtooth  stabilization  or  with  highly  peciked 
heating  profiles.  It  could  also  be  flatter  due  to  naturally  high  elongations,  k  >  2.5. 
Therefore,  a  mean  value  of  ax  =  1.0  appears  consistent  with  the  previous  consider¬ 
ations  and  a  flat  probability  distribution  0.5  <  ax  <  15  is  assumed  as  a  reasonable 
range. 


A. 6  Steady>State  Particle  Balance 

The  steady-state  densities  of  the  various  spedes  present  in  a  burning  plasma 
are  determined  by  a  detailed  particle  beilance  calculation.  In  addition  to  the  con¬ 
sumption  and  production  of  various  spedes  in  a  burning  plasma,  a  realistic  steady- 
state  particle  bcdance  must  consider  the  fuehng  of  the  plasma,  the  presence  of  impu¬ 
rity  concentrations,  the  trapping  fraction  of  the  various  energetic  fusion  products, 
and  the  diffusion  loss  of  charged  particles  out  of  the  plasma.  Initial  assumptions 
must  define  density  and  temperature  parameters,  an  effective  ®He:D  fuel  ratio,  and 
effective  steady-state  impurity  concentrations.  Diffusion  losses  from  the  plasma  are 
based  on  the  calculated  global  energy  confinement  time  rjj,  and  the  required  fuel  in¬ 
jection  rates  are  calculated  to  balance  total  losses  and  maintain  plasma  density.  The 
final  steady-state  particle  balance,  however,  is  constrained  by  the  effective  density 
and  beta  hmits  of  the  plasma  configuration. 

The  steady-state  particle  balance  equations  which  describe  the  diffusion  loss, 
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fuel  injection,  production  and  consumption  of  a  given  species  are  of  the  form 

^  +  '^ninj6ii{<rv)ijJ]f  -  ^n,Ti_f  =  0 

^  ’’’p  ij  * 

where  i  and  j  represent  the  appropriate  species  for  reaction  ij  which  produces  species 

f.  The  term  is  the  diffusion  loss  rate  and  is  a  source  term  representing 

the  fuel  injection  rate  for  species  f.  The  fusion  reaction  rates  {(Tv)ij  which  govern 

particle  production  and  consumption  represent  the  following  five  fusion  reactions 

[7]:  ij=ddl  for  D(D,P)T,  dd2  for  D(D,n)3He,  dhe3  for  D(3He,P)o,  dt  for  D(T,n)a, 

and  tt  for  T(T,2n)Q.  For  like  particles,  ~  1/2  when  i  =  j,  otherwise  Sij  ~  1 

when  i  ^  j-  For  D-T  plasmas  operating  at  10-20  keV,  the  ddl  and  dd2  reactions 

occur  with  nearly  equal  probability  at  a  50:50  ratio.  For  D-^He  plasmas  operating 

at  50-80  keV,  however,  the  ratio  53:47  is  a  more  accurate  depiction  of  the  ddl  and 

dd2  reaction  rates,  respectively. 

The  fractions  of  particles  trapped  on  their  first  orbit  are  modelled  by  adjustable 
parameters  rif,  where  f  represents  3.02  and  14.7  MeV  protons  (written  as  subscripts 
P3  and  P14),  1.01  MeV  tritons  (T),  0.82  MeV  ^He  particles  (®He)  and  the  3.5  and  3.7 
MeV  alpha  particles  (a)  produced  by  the  fusion  reactions.  The  trapping  fractions  of 
fusion  particle  energy  in  the  plasma  depend  on  the  magnit.  de  of  the  plasma  current 
[28].  Since  the  plasma  current  is  generally  much  larger  than  14/ A  MA,  most  of  the 
14.7  MeV  protons  can  be  confined  and  all  of  the  slower  particles  are  assumed  to  be 
confined.  The  default  trapping  fractions  T]f  are  assumed  to  be  1,  except  rjpu  =0.9 
for  the  14.7  MeV  protons. 

All  particle  confinement  times  are  assumed  to  have  the  same  value  of  r^.  If  the 
energy  confinement  time  te  includes  both  convection  and  conduction  power  losses, 
then 

J_  _  2  ^  C 

where  is  the  confinement  time  associated  with  conductive  power  losses  [19].  If 
we  assume  C=l,  a  default  value  that  determines  the  power  loss  associated  with  the 
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outward  flow  of  particles  from  the  core  plasma,  and  a  ratio  =  Tp/xc  =  2,  which 
is  at  the  low  end  of  the  experimental  and  theoretical  indications,  then  the  above 
equation  may  be  rewritten  as 

XjE  Tp  Tp  Tp 

or  simply  Xp/r^  ~  3,  where  the  factor  3  may  be  varied  as  an  input  parameter.  It 
is  important  to  note  here  that  this  is  a  very  simple  approximation,  and  that  global 
parameters  such  as  the  energy  and  particle  confinement  times  do  not  adequately 
describe  the  physics  of  the  D-^He  and  D-D-D  burning  plasmas  [40]. 

The  default  value  for  the  effective  fuel  ratio  /ajp  equal  to  ns/n^j  is  1.0,  which 
assumes  ud  =  713.  An  important  feature  of  D-^He  fuel  cycles  is  that  the  maximum 
value  of  the  fusion  power  density  achievable  is  virtually  independent  of  allow¬ 
ing  optimization  with  respect  to  other  design  parameters  [4].  A  deuterium  lean 
fuel  mixture  has  the  advantage  of  reducing  the  dd  reactions  and  the  production 
of  tritium,  thus  reducing  the  tritium  handling  problem  as  well  as  producing  fewer 
neutrons.  A  deuterium  rich  fuel  mixture  has  the  advantage  of  reducing  the  relative 
electron  density  and  radiation  losses.  An  optimum  fuel  mixture  for  tokamaks  is 
estimated  to  be  about  35%  ^He  and  65%  D  [5].  A  reasonable  range  for  the  fuel 
ratio  is  0.5  <  /axj  <  2.0,  which  represents  a  realistic  operating  domain  within  D-^He 
ignition  space  [4]. 

A  D-^He  plasma  is  assumed  to  contain  reasonable  concentrations  of  represen¬ 
tative  low-Z  (carbon,  oxygen)  and  small  amounts  of  high-Z  (iron,  tungsten,  etc.) 
impurities.  A  low  oxygen  concentration  is  expected  because  tokamaks  with  long 
periods  of  operation  and  hot  walls  have  less  oxygen  than  those  with  cold  walls  or 
low  power  cases,  and  carbon  concentrations  drop  at  high  densities  (ric  ~  l/w^).  The 
iron  concentration  is  typically  very  low  because  little  or  no  Fe-bearing  materials  are 
expected  to  be  used  as  a  first  wall,  but  may  be  considered  representative  of  high- 
Z  impurities  which  may  be  present.  Experimental  observations  give  the  following 
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simple  fits  [35],  normalized  at  (ng)  =  0.7  x  which  are  used  as  defaults: 

nc/n,  =  [0.9  +  0.6(0.7e20/nef  •®]% 


=  01% 

™Fe/n«  =  [0.05(0. 7e20/ne)^®]% 

Seven  simultaneous  equations  may  now  be  written  which  will  solve  for  seven 
unknowns:  rip,  no,  nj,  na,  tIq,  Sd,  and  S3.  Particle  balance  equations  for  the 
five  species,  a  charge  balance  equation  (which  includes  a  term  for  the  neutral  beam 
density,  and  the  fuel  ratio  equation  are  required  for  solution  of  the  steady- 

state  particle  balance.  These  equations  are: 


rip 


—  —  ^A'lnjj{(Tv)^irfp3  +  n£,n3{(rv)a},g3r)pi4 


np) 


—  So  —  0.47njp(<n;)j^i  —  0.53np(oT;)^2  ~  ~  ^d'^t{'^)cU 


rif 


=  St  +  O.A7np{(rv)jair}T  -nonT{(rv)^  -np{<rv)tt 


71.3 


-  —  ^3  0.53n£,(<rv)tW2J7’He  ~  ”D^3(<^)(i/ie3 


—  —  rion3{<rv)j}te3Va  +  nonT{<^)dt<la  +  np{<rv)tt 

’’’p 

rig  =  {np  -f  np  -H  ny  -f  nj,eom)  +  2(713  +  ^a)  +  +  StIq  +  267ljre 


The  fusion  reaction  rates  {<rv)ij  (where  cr  is  the  reaction  cross  section,  v  is 
the  relative  speed  of  the  reactants,  and  the  brackets  indicate  averaging  over  the 
Maxwelhan  velocity  distribution  of  the  reacting  species)  are  volume  averaged  using 
the  density  and  temperature  profiles  where  =  {rii),  rij  =  {rij),  and 

(<n>)„(a„,a,,  (T))  =  2(1  +a„f  [' p(l  -  +  ar)(l  -  P^r^]dp 

Jo 

where  {ov)ij  is  the  temperature  dependent  Maxwellian  reaction  rate  [41]. 
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Disregarding  the  second  and  fourth  equations  for  now,  the  other  five  equations 
are  used  to  solve  for  np,nT,nz  and  in  terms  oinj).  Assuming  tritium  recycling, 
in  which  tritium  diffusion  losses  are  injected  back  into  the  pl2isma  at  an  equal  rate, 
allows  for  the  substitution  St  =  nxlrp  in  the  third  equation.  No  additional  tritium 
fueling  of  the  plasma  is  assumed.  Recycling  of  tritium  hardens  the  neutron  spectrum 
^d  does  not  do  much  to  improve  the  attainable  power  density,  but  does  minimize 
the  tritium  storage  requirement  (with  its  inherent  radioactivity  and  weapons  pro¬ 
liferation  risks).  Letting  ns  =  nufsD  and  n*  —  ntgam  ~  6nc  —  8no  —  26n/’e  =  C,  a 
constant,  the  original  seven  equations  are  reduced  to  the  following  four  equations 
with  four  unknowns: 

—  =  0A7nj^{ffv)jaiVP2  +  ^DfaD{o^)dht3!flPu 
0  =  0.47n^(<n;)^ii7r  —n^nT^ov)^  -n^{(rv)u 
^  =  nl,f3D{ov)dheS‘na  +  nDnT{ov)jtVa  +  nl{<rv)uVa 

~p 

C  —  np  —  (1  •+•  2/3ij)np  —  nx  —  2na  =  0 

Now  np,  np  and  n,,  are  calculated  in  terms  of  nu  only,  with  all  other  terms 
representing  constants. 


np  =  rp(0.47(<n;)rfrfiqp3  -f /3D(<n;)jhe3’7i»i4)np 


np  = 


+  l-SBn]j{<rv)ddi{<ro)uVT 
2{<rv)tt 

na  =  rp{nljfiD{<rv)dhe3Ve,+nDnT{<ro)dtVa+np{av)ttVa) 


Note  now 

C  —  np  —  (1  -f  2/3x>)ni;  —  np  —  2no  =  0 

is  a  function  of  n©  only  and  may  be  solved  numerically  for  the  steady-state  deuterium 
density  in  the  plasma.  STORAC  solves  for  the  zero  of  the  function  f(n£))=0  and 
returns  the  value  for  n^j.  The  other  steady-state  particle  densities  are  then  calculated 
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using  the  previously  derived  equations.  The  second  and  fourth  equations  of  the 
original  set  are  now  used  to  calculate  the  source  terms,  Sd  and  Sz,  representing  the 
fuel  injection  rates  required  to  maintain  the  steady-state  particle  balance. 

A.7  Density  Limit 

The  density  hmit  imposes  an  upper  limit  on  the  plasma  edge  density.  A  ’’weak” 
hmit  is  usually  assumed,  implying  that  no  a  priori  density  hmit  is  to  be  imposed.  An 
optimistic  estimate  of  what  is  achievable  under  best  conditions  is  probably  provided 
by  the  boundary  of  the  accessible  edge  physics  parameter  space.  Tentatively,  the 
Borrass  density  hmit  for  an  ITER  hke  device  given  by  [35]  is 

where  is  the  critical  plasma  density  at  the  separatrix,  =  q^(95%), 

P ±{MWfm^)  is  the  mean  heat  flux  across  the  separatrix,  and  Cwl.8  for  ITER-hke 
conditions.  STORAC  uses  a  modified  Borrass  model  for  the  spherical  tolmmak  of 
the  form 

^  0  1425  X 

We96 

where  the  factor  0.1425  x  10^  accounts  for  the  density  in  the  scrape-off  layer  relative 
to  the  average  plasma  density.  The  density  hmit  may  be  exceeded,  which  has  been 
observed  experimentally  in  discharges  with  pellet  injection  or  strong  additionad  heat- 
ing  [36].  To  allow  for  shghtly  higher  densities  during  H-mode  operation,  the  density 
hmit  constreiint  is  modified  by  y/2  [2?]  and  imposed  as  a  default. 

A. 8  Beta  Limit 

The  plasma  beta,  is  defined  as  the  ratio  of  kinetic  pressure  to  confining 
magnetic  pressure  and  is  a  measure  of  the  efficiency  with  which  the  magnetic  field 
confines  the  plasma.  The  total  volume  average  plasma  beta  is  the  sum  of  the  thermal 
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plasma  {/3th),  fast  particle  {/Sfast),  and  energetic  beam  (/Sseom)  contributions.  These 
components  are: 


13th  = 


{rieT^  +  njTj) 

By2fio 


_  2nf{Efa,t)/^ 

BV2/io 


Pheam  — 


2n 


beam beam 


/3 


BV2/io 

where  B  is  the  total  magnetic  field  in  the  plasma  approximated  as  (B  j  +  Biyi\  ,.0 
is  the  permeability  of  free  space  (47r  x  10"^  H/m),  {Ejatt)  is  the  average  fast  particle 
energy  approximated  as  O.SOEfatt  (about  1/3  of  the  fast  particle’s  birth  energy), 
and  the  beam  density  and  energy  are  defined  by  current  drive  requirements. 

The  plasma  beta  can  be  limited  by  either  the  Troyon  beta  hmit  (/3y,  appro¬ 
priate  for  the  first  stability  regime)  or  the  e/0p  limit  (/S,,  appropriate  for  the  second 
stabihty  regime),  where 

IpiMA) 


/9/<10-"xp/ 


and 


a{m)Bt{T) 
5 


For  steady-state  operation,  a  Troyon  co^dent  3.0  <  gf  <  3.3  seems  to  represent 
the  most  likely  region  of  operation,  although  the  DIII-D  tokamak  has  operated 
stably  with  gj  as  high  as  3.5  [39].  The  default  g,  value  for  the  maximum  e/3p  is  0.6 
for  operation  in  the  1st  stability  regime  and  1.0  for  operation  in  the  second  stability 
regime,  e  is  the  inverse  ^lspect  ratio  and  the  poloidal  beta  is 


a  —  _  a  (  B  \ 

5^V2/io  ~^[bJ 


Due  to  plasma  paramagnetism  effects  and  the  definition  of  the  total  volume 
averaged  plasma  beta  used  in  this  model,  beta  values  may  exceed  unity  as  ej3p  values 
approach  or  exceed  unity  [27].  The  e/3p  limit  is  considered  to  be  a  more  important 
constraint  than  the  Troyon  beta  hmit,  but  there  is  much  speculation  regarding  how 
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high  e/9p  can  safely  go.  However,  requiring  the  bootstrap  current  to  be  less  than  the 
total  plasma  current  for  equilibrium  (see  A.4)  may  represent  a  t3rpe  of  beta  hmit. 
This  was  not  included  in  any  calculations  for  this  study,  but  could  be  considered  in 
future  studies. 


A.9  Global  Plasma  Power  Balance  Equations 

The  general  global  plasma  power  balance  equation  for  D-^He  plasmas  with 
electron  temperature  Tg  and  ion  temperature  Ti  =  Ti,  =  given  by  [28]  is 

^  +  n-a  +  ngl'n)KTi  =  Pf  +  Poh  +  Poux  -  Peon  -  Prod  -  P« 

where  7^  =  TJTg  is  the  ion  to  electron  temperature  ratio,  K  is  the  Boltzmann 
constant,  Pp  is  the  fusion  power  density,  Poh  is  the  ohnoic  heating  power  density, 
Paux  is  the  auxiliary  heating  power  density,  Pam  is  the  conduction  and  convection 
loss  term,  P^ad  is  the  bremsstrahlung  and  synchrotron  radiation  loss  term,  and  Pje 
is  the  equilibration  power  between  ions  etnd  electrons  when  Tj  and  Tg  are  not  equal. 
This  study  assumes  Ti  —Tg,  so  Pj*  =  0.  All  background  ions  are  assumed  to  be  at 
the  same  local  temperature,  Ti{p\  and  to  have  a  common  profile.  In  steady-state, 
the  term  equals  zero,  thus 

Pf  +  Poh  +  Paux  =  Peon  +  Prod 

The  global  plasma  power  balance  constraint  equation  used  in  STORAC,  which 
requires  an  ignition  margin  to  equal  total  power  sources  divided  by  total  power  losses, 


'  Peon  Prod-c 
\Pf  -I-  Poh  +  Pot 


X  IM  <  e 


where  c  ~  0  is  an  error  tolerance,  usually  10 or  10“^.  IM  is  the  ignition  margin, 
assumed  to  be  1  by  default.  Only  the  core  radiation  power  loss  (Prod-c)  is  included 
in  the  global  power  balance  equation.  The  edge  radiation  is  assumed  to  be  lost  to 


the  divertor  region. 
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The  fusion  power  density  (in  MWfm}),  calculated  by  summing  over  all  fusion 
reactions,  is; 

Pp  =  [0.47n|j(<n;)d(fi(1.0l7/j’ +  3. 02^jc>3)  +  0.53n^(<rv  )dj2(0. 827/3^^5 )  + 

nonz{<rv)jitei{Z.lr}a  +  14.777pi4)  +  nDn7'(<n;)eit(3.57/a)  + 

np{<Tv)ti{1.26irja)]  x  1.6  x  10“^® 

The  ohmic  heating  power  density  is: 

where 

Coh  «  9.33  X  10-^Z,ff[{l  +  a„)/(l  +  a„  +  ar)]"/*(l  +  3ar/2) 

and 

"yNc  —  4*3  ~  0.6i4 

is  a  neoclassical  resistivity  enheincement  factor  (ITER  definition)  [35].  Since  ohnoic 
heating  at  thermonucleu  temperatures  is  a  very  slow  process  and  there  is  no  loop 
voltage  for  steady-state  operation,  P^h  for  steady-state  operation  is  negligible. 

The  auxileiry  power  density  is  a  function  of  the  injected  power  to  the  ions  and 
electrons  required  for  current  drive  and  plasma  heating.  Neutred  beam  injection  is 
used  as  the  default  method  for  current  drive,  so  Pa„  is  a  function  of  the  neutral 
beeim  power  required  for  noninductive  current  drive.  The  current  drive  subroutine 
in  STORAC  determines  the  neutral  beam  current  drive  efficiency,  coiilomb  loga¬ 
rithm  for  ion-ion  collisions,  and  the  stopping  cross-section  for  a  hydrogen  beam  in  a 
fusion  plasma.  The  original  equations  used  in  STORAC  calculate  the  neutral  beam 
efficiency  using  modified  fitting  functions  from  the  ITER  neutral  beam  model  [35]. 
These  were  further  modified  for  a  D-^He  plasma  background  by  replacing  tritium 
with  ®He  where  appropriate.  The  original  Tio  scalings  (Tio  =  {T]/10keV)  were  re¬ 
tained,  which  may  not  be  appropriate  for  the  higher  operating  temperatures  of  a 
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D-^He  plasma.  Although  the  neutral  beam  model  is  somewhat  speculative,  it  is  still 
useful  for  determining  auxihary  power  requirements  and  their  estimated  costs. 

The  conductive  and  convective  losses  (transport  losses)  are  calcidated  using 
global  energy  confinement  scalings.  These  scalings  are  usually  empirical,  and  predict 
the  global  energy  confinement  time  (te)  as  a  function  of  device  parameters.  The 
power  density  of  the  transport  losses  is 

p_  =  =  2.404  X  l0-”((n,r,}  +  (n.-T,))  ^ 

Te  Te 

where  Wth  is  the  plasma  thermal  energy  content.  For  a  pure  homogeneous  pleisma, 
this  is  simply  =  2nKT  xVol. 

The  1990  ITER  H-mode  power  scaling  was  assumed  for  this  study.  It  is 

where  P^t  is  the  total  power  injected  into  the  plasma  from  outside  including  all 
heating  and  current  drive  methods  and  Ai  is  the  average  atoiriic  mass  number  of 
the  plasma  ions.  The  parameter  fx,,  or  H-&.ctor,  is  a  confinement  time  enhance¬ 
ment  factor  relative  to  L-mode.  Improved  confinement  is  observed  particularly  in 
machines  with  a  divertor  when  the  power  input  exceeds  a  certain  threshold  value, 
with  1  <  <  4  representing  a  reasonable  range  of  values.  At  the  end  of  each  run, 
STORAC  calculates  fi  required  to  maintain  power  bzilance. 

The  radiation  power  term  includes  bremsstrahlung  and  synchrotron  radiation 
calculated  for  the  core  and  edge  plasma.  Only  the  core  radiation  power  loss  is 
included  in  the  plasma  power  balance.  The  edge  radiation  is  subtracted  from  the 
power  flow  to  the  divertor  region.  Bremsstrahlung  (or  braking)  radiation  in  fusion 
plasmas  is  due  to  collisions  between  a  Maxwellieui  sea  of  electrons  with  ions  of  charge 
Zi-  The  electron-electron  and  electron-ion  bremsstrahlung  loc2il  power  density  given 
by  [7]  is 

=  5.355  X  1.55  x  10'^,  +  7.15  x  lO"®?’/ 
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which  incorporates  first-order  e-i  relativistic  and  spin  corrections,  e-e  emission  (non- 
relativistic),  and  a  Bom  correction  [7].  Z^ff  is  the  effective  charge  of  the  plasma  ions, 
where  Ze//  =  The  edge  radiation  power  includes  separate  calculations 

for  specified  impurity  density  fractions  and  an  assumption  of  coronal  equilibrium. 
Synchrotron  (or  cyclotron)  radiation  is  due  to  the  acceleration  of  a  charged  particle 
moving  perpendicular  to  a  magnetic  field  as  it  executes  a  gyro-orbit,  which  results  in 
emission  of  electromagnetic  radiation.  The  synchrotron  global  power  density  given 


by  [7]  is 


=  6.214  X 


where  the  assumed  Tmbnikov  absorption  correction  is 


=  7.78  X  10-» 


The  wall  reflectivity  is  assumed  to  be  0.9,  accounting  for  wall  deterioration  and 
penetrations  of  a  smooth  metallic  wall.  Corrections  for  plasma  diamagnetism  are 
approximated  by  substituting  Bo(l  —  for  the  vacuum  field  Bq.  Technically, 
this  diamagnetic  correction  is  inappropriate  when  the  plasma  displays  paramag¬ 
netic  behavior,  but  synchrotron  radiation  is  expected  to  be  neghgible  at  high  beta 
regardless.  The  total  radiation  power  is  then  the  sum  of  the  bremsstrahlung  and 
synchrotron  radiation  terms. 

K  a  separate  power  balance  for  ions  and  electrons  is  used,  then  the  equilibration 


power  between  ions  and  electrons  is 


Pie  =  2.4  X  10  *^CiJ,nkien\Zef] 


{Ti-Te) 


n  ^  _ (1  +Qn)^ _ 

*'  (2a„  -  0.5ar  +  1)(1 -h  Or)'/' 


where 
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is  a  volume  averaged  profile  factor.  Of  course,  if  T,-  =  T^,  then  Pie  =  0.  An 
increase  in  the  ratio  of  Ti/Tg  wiU  lead  in  most  cases  to  increased  fusion  reaction 
power  compared  to  the  radiation  power  loss  [7].  This  potential  enhancement  of 
plasma  performance,  and  the  physics  factors  affecting  the  temperature  ratio,  would 
be  worth  investigating  in  future  studies. 


APPENDIX  B 

Reference  Cases  ST-lB  and  ST-2B 


The  reference  cases  used  for  this  study  have  been  designated  ST-IB  auttd  ST-2B, 
respectively.  Both  cases  assume  an  edge  safety  factor  ~  ^)  '^hich  is  a  compromise 
between  the  physics  advantages  of  high-q  operation  and  the  engineering  advantages 
of  low-q  operation,  as  discussed  in  Chapter  3.  The  parametric  investigation  of  e/dp 
values  established  a  linkage  between  the  1st  and  2nd  stability  regimes,  allowing 
further  parametric  investigations  to  be  made  based  on  just  the  ST-2B  reference 
case.  The  anticipated  effects  of  specific  parameter  variations  on  reactor  performance 
relative  to  either  reference  case  are  expected  to  be  comparable. 

Below  is  a  brief  review  of  the  assumptions  adopted  for  the  ST- IB  and  ST-2B 
reference  cases,  as  discussed  in  Chapter  2: 

•  Aspect  ratio,  Rofa  >  1.2 

•  Plasma  elongation,  =  2. 1(1.0  -f-  0.44c*-^) 

•  Plasma  triangularity,  8x  =  0.53(1.0  ■¥  0.77c®  ®) 

•  e/3p  =  0.6  for  ST-IB,  1.0  for  ST-2B 

•  Central  safety  factor  q(0)=1.0  for  ST-IB,  2.0  for  ST-2B 

•  Temperature  profile  factor,  ax  =  1-0 

•  Density  profile  factor,  a„  =  0.5 

•  Effective  fuel  ratio,  ns/njj  =  1.0 

•  Confinement  time  ratio,  Tp/r^  =  3.0 

•  1990  ITER  H-mode  power  scaling  for 

•  Bootstrap  current  fraction,  ht/Ip  —  0-90 

■  Centerpost  current  density,  J^p  <  75MAfm^ 

•  Thermal  conversion  efficiency,  rith  =  0.40 

•  Direct  energy  conversion  efficiency,  i^jec  =  0.50 

•  Net  electric  power  output,  1000  MWe 

•  Plant  operating  life,  30  years 

•  Plant  capacity  factor,  0.75 

•  Level  of  safety  assurance,  LSA  =  2 


94 


95 


Parameter 

ST-IB 

ST-2B 

Error  of  Constraint  Equations: 

Beta 

-1.58e-07 

2.99e-06 

Global  power  bedance 

-2.02e-07 

3.79e-06 

Beam  ion  density 

1.84e-07 

-3.92e-06 

Field  at  coil 

-4.06e-12 

8.83e-09 

Radial  build 

-7.11e-15 

-2.83e-09 

Centerpost  temperature 

2.00e-07 

-1.32e-06 

Particle  confinement  time 

-2.47e-06 

4.65e-05 

Density  limit 

-2.62e-07 

4.40e-06 

Troyon  beta  limit 

4.26e-14 

1.42e-14 

limit 

-8.53e-14 

1.85e-13 

Beam  energy 

-3.61e-09 

-1.04e-08 

Net  electric  power 

2.06e-07 

-2.94e-06 

Average  error  of  constraints 

2.52e-06 

4.73e-05 

Plasma  Geometry  Parameters: 

Major  radius,  Rq  (m) 

2.154 

1.912 

Minor  radius,  a  (m) 

1.795 

1.594 

Aspect  ratio,  A  =  Ro/a 

1.200 

1.200 

Elongation,  X-point 

2.665 

2.665 

Elongation,  95%  surface  K95 

2.386 

2.386 

Triangularity,  X-point  Sx 

0.766 

0.766 

Triangularity,  95%  surface  Sgs 

0.501 

0.501 

Plasma  surface  area,  5p(m*) 

247.2 

194.9 

Plasma  volume,  Vp(m^) 

297.4 

208.2 

Current  and  Field  Parameters: 

Plasma  current.  Ip  (MA) 

112.404 

84.185 

On-axis  toroidal  magnetic  field,  Bq  (T) 

1.955 

1.649 

Average  poloidal  magnetic  field.  Bp  (T) 

6.487 

5.472 

Total  magnetic  field,  B  (T) 

6.775 

5.715 

On-axis  safety  factor,  qo 

1.000 

2.000 

Edge  safety  factor, 

5.000 

5.000 

95%  surface  safety  factor,  995 

2.218 

2.218 

Parameter 


ST-IB 


ST-2B 


Beta  Parameters: 

Pleisma  beta,  ^ 

0.660 

1.100 

Poloidal  beta, 

0.720 

1.200 

Fast  particle  beta,  ^fatt 

9.52e-02 

1.74e-01 

Neutral  beam  beta, 

6.32e-03 

8.12e-03 

Troyon  coefficient,  gj 

2.061 

3.435 

0.600 

1.000 

Beta  limit,  /3 j  and 

1.057 

1.100 

Temperature  and  Density  Parameters: 

Temperature  profile  factor,  aj 

1.000 

1.000 

Density  profile  factor,  a„ 

0.500 

0.500 

Electron  temperature,  (Tg)  (keV) 

59.877 

62.937 

Ion  temperature,  (Tj)  (keV) 

59.877 

62.937 

Electron  n-temperature,  [Tg]  (keV) 

71.852 

75.524 

Electron  density,  (ng)(10“Tn“®) 

5.086 

5.640 

Ion  density,  (ni)(10“m~^) 

3.776 

4.223 

Deuterium  density,  (n£))(10”m"®) 

1.526 

1.692 

Helium-3  density,  (n3)(10*°m“^) 

1.526 

1.692 

Tritium  density,  (nr)(10^®m~®) 

1.433 

1.682 

High-Z  impurity  density,  {nz)(10^®Tn~®) 

5.106 

5.657 

Alpha  particle  density,  (nQ)(10^°m~®) 

0.344 

0.401 

Proton  density,  {np)(10*°m“®) 

0.312 

0.363 

Effective  charge,  Z^ff 

2.228 

2.244 

Mass- weighted  Zg/y 

0.711 

0.718 

Deuterium  fuel  rate,  5ij(10*’s“^Tn“®) 

0.183 

0.229 

Helium-3  fuel  rate,  53(10’®^ 

0.174 

0.217 

Fusion  Power  Parameters: 

Charged  fusion  power,  Pp  (MW) 

2585 

2348 

Beam  fusion  power,  Pp-Beam  (MW) 

8.970 

6.196 

Neutron  fusion  power,  (MW) 

165.4 

150.1 

Neutron  wall  load  {MWjm}) 

0.535 

0.616 

Fraction  of  power  to  electrons 

0.306 

0.287 

Fraction  of  power  to  ions 

0.694 

0.713 
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Parameter 

ST-IB 

ST-2B 

Plasma  Power  Balance  Parameters: 

Ohmic  heating  power,  Poh  (W) 

0.881 

0.520 

Bremsstreihlung  radiation  power,  Pj,,  (MW) 

1329 

1191 

Synchrotron  radiation  power,  (MW) 

4.684 

.0001 

Wall  reflectivity, 

0.900 

0.900 

Edge  radiation  power,  Pune  (MW) 

19.69 

39.96 

Electron  transport  power,  Pt,_e  (MW) 

748.8 

682.1 

Ion  transport  power,  P*,  j  (MW) 

556.0 

510.8 

Injected  power  to  electrons,  Pinj-e  (MW) 

26.72 

17.61 

Injected  power  to  ions,  Pmj-i  (MW) 

26.12 

18.48 

Power  to  divertor,  P*„ 

1255 

1153 

Divertor  heat  load  (MW/m^) 

1.880 

2.180 

Divertor  Hfe  (years) 

2.120 

1.830 

ITER-H90P  confinement-time  multiplier,  r 

2.582 

3.633 

Globed  confinement  time,  (s) 

3.487 

3.124 

Particle  confinement  time,  Tp  (s) 

10.46 

9.371 

Lawson  parameter,  nrf;(10*^s/Tn®) 

1.774 

1.762 

NBI  Current-Drive  Parameters: 

Steady-state  power  required  (MW) 

52.84 

36.09 

Heating  power  (MW) 

3.375 

1.002 

Plasma  Q- value,  Qp  =  PfIPcd 

48.92 

65.06 

Bootstrap  current  fraction,  htllp 

0.900 

0.900 

Non-inductive  current  fraction 

0.100 

0.100 

Neutral  beam  efiidency,  i^seom  (A./W) 

0.227 

0.240 

Neutral  beam  gamma,  AfW  —  m*) 

2.489 

2.587 

Neutral  beam  power,  Pbeam  (MW) 

52.85 

36.09 

Neutrad  beam  energy,  Ebeam  (MeV) 

5.026 

4.995 

Neutral  beam  current,  /be«m  (A.) 

10.51 

7.226 

Beam  energy  fraction  to  ions 

0.494 

0.512 

Neutral  beam  shine  through 

2.48e-3 

2.48e-3 

Beam  decay  lengths  to  Rq 

3.000 

3.000 
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Parameter 

ST-IB 

ST-2B 

Radial  Build  Paxameters  (m): 

Device  centerline 

0.000 

0.000 

Centerpost 

0.299 

0.259 

First  wall  -  inboard 

0.319 

0.279 

Scrape  off  layer  -  inboaxd 

0.359 

0.319 

Plasma  centerline,  Rq 

2.154 

1.912 

Plasma  edge 

3.948 

3.506 

Scrape  off  layer  -  outboard 

3.988 

3.546 

First  wall  -  outboard 

4.008 

3.566 

Shield  -  outboard 

4.508 

4.066 

Gap 

5.043 

4.481 

TF  coil 

5.342 

4.740 

Vertical  Build  Parameters  (m): 

Plasma  midplane 

0.000 

0.000 

Plasma  top 

4.783 

4.247 

Scrape  off  layer 

7.924 

7.036 

Divertor  structure 

8.024 

7.136 

Shield 

8.524 

7.636 

Gap 

8.687 

7.799 

TF  coil 

8.986 

8.058 

Port  requirement  for  beams  (m) 

0.574 

0.534 

Port  width  (m) 

1.592 

1.439 

Centerpost  Parameters: 

Centerpost  current  density  (MA/m^) 

75.00 

75.00 

Average  resistivity  (10“*  Ohm-m) 

2.364 

2.358 

Resistive  power  loss  (MW) 

230.2 

154.1 

Nuclear  heating  (MW) 

74.51 

64.73 

Centerpost  coolant  fraction 

0.106 

0.105 

Coolant  channel  diameter  (cm) 

1.224 

1.190 

Coolant  channel  length  (m) 

17.37 

15.60 

Coolant  flow  speed  (m/s) 

20.00 

20.00 

Number  of  coolant  tubes 

800.3 

628.8 

Inlet  coolant  temperature  ("C) 

40.00 

40.00 

Paxameter 

ST-IB 

ST-2B 

Centerpost  Parameters,  continued; 

Coolant  temperature  rise  (°C) 

123.9 

120.0 

Average  centerpost  temperature  (“C) 

116.6 

115.2 

Peak  centerpost  temperature  (“C) 

180.1 

176.7 

Coolant  inlet  pressure  (MPa) 

4.801 

4.451 

Coolant  pressure  drop  (MPa) 

3.712 

3.448 

Pump  power  required  (MW) 

8.738 

6.024 

Centerpost  volume  (m^) 

15.38 

10.35 

Weight  of  centerpost  (10®  kg) 

0.122 

0.083 

TF  Coil  Parameters: 

Number  of  TF  coil  legs 

16 

16 

Outer  leg  current  density  (MA/m^) 

4.369 

4.799 

Coil  resistance  (10  “®  Ohm) 

2.391 

3.123 

Resistive  power  loss  (MW) 

66.2 

48.6 

TF  coil  current  (MA) 

21.06 

15.77 

Peak  field  at  TF  coil,  (T) 

14.09 

12.19 

Ripple  at  plasma  edge  (%) 

1.50 

1.50 

Stored  energy  per  coil  (GJ) 

0.136 

0.069 

Vertical  force  -  inner  leg  (MN) 

4.307 

2.438 

Radial  stress  (MPa) 

118.0 

88.25 

Coil  volume  (m®) 

8.679 

5.271 

Weight  of  outer  legs  (10®  kg) 

0.989 

0.601 

TF  Power  Supply  Parameters: 

Bus  current  density  {MAfm^) 

1.250 

1.250 

Maximum  impedance  (Ohm) 

4.207e-05 

5.507e-05 

Bus  resistance  (Ohm) 

3.799e-06 

5.073e-06 

Bus  length  -  all  coils  (m) 

160.0 

160.0 

Resistive  power  loss  (MW) 

6.6 

4.9 

Peak  voltage/coil  (V) 

3.46 

3.39 

Peak  power  (MW) 

30.30 

20.75 

Parameter 

ST-IB 

ST-2B 

- - - - - - -  1 

PF  Coil  Parameters: 

Number  of  PF  coil  groups 

2 

2 

Circuit  1,2  radius  from  Ro  (m) 

0.78 

0.69 

Circuit  1,2  distance  from  midplane  (m) 

±8.02 

±7.14 

Circuit  1,2  field  at  coil  (T) 

27.5 

23.5 

Circuit  1,2  current  (MA) 

45.14 

33.81 

Circuit  1,2  current  density  (Af  A/m*) 

5.624 

0.673 

Circuit  3.4  radius  from  (m) 

5.84 

5.24 

CircTiit  3,4  distance  from  midplane  (m) 

±2.15 

±1.91 

Circuit  3,4  field  at  coil  (T) 

27.5 

23.5 

Circuit  3,4  current  (MA) 

-44.96 

-33.67 

Circuit  3,4  current  density  (Af  A/m*) 

3.949 

3.772 

Bus  length  -  all  circuits  (m) 

314.5 

310.6 

Resistive  power  loss  (MW) 

0.31 

0.31 

Majdmum  PF  coil  voltage  (kV) 

5.0 

5.0 

Outer  Shield  Parameters: 

Shield  thickness  (m) 

0.50 

0.50 

Shield  volume  (m*) 

190.8 

154.9 

Weight  of  shield  (10®  kg) 

1.116 

0.906 

Nuclear  heating  (MW) 

1523 

1366 

Reactor  Power  Parameters: 

Charged  fusion  power  (MW) 

2585 

2348 

Neutron  fusion  power  (MW) 

165 

150 

Neutron  power  multipher 

1.270 

1.270 

Injector  wall  plug  power  (MW) 

105.7 

72.2 

Vacuum  heat  (MW) 

0.50 

0.50 

Tritium  process  heat  (MW) 

1.00 

1.00 

TF  coil  resistive  power  (MW) 

303.0 

207.5 

Coolant  pump  power  (MW) 

8.7 

6.0 

Power  loss  to  holes  (MW) 

224.8 

201.2 

Primary  heat  (MW) 

2151 

1942 

Secondary  heat  (MW) 

59.1 

45.2 

Parameter 

ST-IB 

ST-2B 

Reactor  Parameters,  continued: 

Thermal  conversion  efficiency  (%) 

0.40 

0.40 

DEC  conversion  efficiency  (%) 

0.50 

0.50 

Gross  electric  power  (MWe) 

1488 

1353 

Net  electric  power  output  (MWe) 

1000 

1000 

BOP  recirculating  power  fraction 

0.047 

0.050 

Total  recirculating  power  fraction 

0.330 

0.260 

Cost  of  electricity  (mill/kWh) 

88.90 

73.68 

Cost  Accounts  Output  ($M): 

211  Site  and  facilities 

10.00 

10.00 

212  Reactor  building 

66.84 

59.88 

213  Other  buildings 

90.78 

85.49 

21  Total 

167.62 

155.37 

2211  First  wall 

33.31 

26.41 

2213  Shield 

59.16 

48.04 

2215  Divertor/DEC 

48.57 

38.30 

221  Total 

420.27 

238.82 

2221  TF  magnet  assemblies 

112.90 

69.76 

2222  PF  magnet  assembhes 

285.13 

160.81 

222  Total 

398.02 

230.57 

223  Total,  current  drive  system 

303.83 

220.07 

224  Total,  vacuum  system 

73.94 

69.42 

2251  TF  power  supphes 

178.02 

134.12 

2252  PF  power  supphes 

32.33 

25.93 

225  Total 

210.35 

160.05 

226  Total,  heat  tramsport  system 

199.04 

180.91 

227  Totad,  fuel  handling  system 

259.22 

251.43 

228  Total,  instrumentation  and  control 

182.88 

177.38 

229  Total,  maintenance  equipment 

151.35 

146.80 

22  Total 

2198.92 

1675.43 

23  Turbine  plant  equipment 

391.31 

359.60 

24  Electric  plant  equipment 

116.94 

109.44 

25  Miscellaneous  plant  equipment 

61.96 

57.98 

26  Heat  rejection  system 

84.16 

76.01 

Parameter 


ST-IB 


ST-2B 


Cost  Accounts  Output,  continued  ($M); 

2  Total  plant  direct  cost 

3020.89 

2433.84 

Indirect  costs 

847.66 

682.94 

Contingency  costs 

669.26 

539.20 

Interest  during  construction 

788.56 

635.32 

Total  capital  cost 

5326.37 

4291.30 

Miscellaneous  Parameters: 

Intercoil  structure  mass  (10®  kg) 

1.298 

0.650 

PF  coil  support  mass  (10®  kg) 

7.679 

3.396 

Gener2il  support  mass  (10®  kg) 

0.030 

0.021 

Total  structure  mass  (10®  kg) 

9.007 

4.067 

Reactor  building  volume  (lO^m®) 

73.98 

61.97 

APPENDIX  C 

Edge  Safety  Factor  Scan  of  ST-2A 

This  appendix  contains  all  the  output  plots  generated  by  STORAC  for  the  edge 
safety  fector  scan  of  the  ST-2A  2nd  stability  reference  case  as  was  varied  from 
7.68  to  3.5.  These  are  representative  of  the  plots  generated  from  the  plot  files  of 
every  STORAC  run  which  uses  the  scan  option.  The  scan  routine  stores  in  an  array 
many  of  the  output  quantities  calculated  in  the  optimization  process,  and  prints 
them  into  a  separate  file  which  may  be  read  by  a  plotting  routine. 

Systems  codes  are  often  used  to  perform  scans  in  order  to  determine  the  sen¬ 
sitivity  of  results  to  different  input  assumptions.  STORAC  is  capable  of  performing 
up  to  10  successive  calculations  for  a  prescribed  variable,  where  each  calculation 
iterates  from  the  previous  solution.  The  "smoothness”  of  a  given  plot  depends  on 
the  residual  error  of  the  previous  solution  and  the  step-size  of  the  scan  variable.  If 
either  the  residual  error  or  the  step-size  is  too  large,  successive  iterations  can  pro¬ 
duce  poor  solutions  or  fail  to  execute.  If  the  poor  solution  has  a  residual  error  in 
the  range  of  10  —  10"®  ,  then  successive  iterations  could  still  produce  successful 

optimizations  rather  than  get  worse.  For  example,  the  "dent”  in  the  curve  in  Figure 
3.18  at  Jep  =  'IbMAIrn}  (ironically,  the  assumed  value  for  the  reference  cases)  shows 
how  the  scan  can  recover  from  a  bad  iteration  and  produce  an  otherwise  smooth 
curve. 


103 


Toroidal  Field  On>Axia,  £o  ('I')  Major  Iladiua,  Jto  (m) 


2.50 

2.25 
2.00 
1.75 

1.50 

1.25 

1,00 


3  5 


Edge  Sa&tT'  Factor, 


Edge  Sa&ty  Factor, 


1.2Q 

1.  15 

1.  10 

--  1.05 

ts 

(S 

I  1.00 

s 

^  0.95 

Edge  Saiet]r  factor, 


•• 

O 

X 


o 


*s 


Edge  Safety  Factor, 


Plaima  Current,  //,  (MA)  Uoir=*  Ilpwr=» 


5 


1.00 


0.95  - 

a. 

*^Q.9Q  -  . 

m 

& 

2-0.85  - 

"o 

m  0.80  - ' - ^ ^ 

3  5  7 

Edge  Sail^  Factor, 


Edge  Safsty  Factor, 


Total  Oapital  Goal  ($11)  Temperature,  (T.)  (xlO  keV) 


8.25 

8.00 

7.75 
7.50 

7.25 
7.00 

6.75 


5.5 
5.0 

4.5 
4.0 

3.5 


3  5  7 


Edge  Sa&ty  Factor,  q^ 


Dig  Q,  PfjPao  Direci  Oo*l  (|B) 


3.25 
3.00 
2.75 
2.50 

2.25 

2.00 

3  5  7 

E«lge  Sa£Bt7  Factor, 

67.5 

66.5 

65.5 

64.5 

Edge  Safety  Factor,  q^ 


3  5 

Edge  Safirtj  Fictor,  ^ 


OooUni  Channel  Radius  (xlO~*in)  Ream  Energy,  Eww*  (MeV) 


108 


Edge  Sa&ty  Factor,  9^ 


Edge  Sa&t7  Factor,  9^ 


20 . 10 


Edge  Sattt7  Factor,  ^ 


Cuneni  DeniUy,  Jcf  (MA/m*) 


3  5 

Edge  Sa£rt7'  Factor,  g* 


.  10 
75.05 
75.00 
74.95 
74.30 


Edge  Sa£st7  Factor, 


3: 

S 


a. 

E 


a. 


Edge  Safety  Factor, 


'Tl.OO  10 
?1.00  10 
r 1 . 0  0  0  0 

41 

*^1.0000 

•S 

ll.OOOO 
*0 . 9938 


3  5 


Edge  Safety  Factor,  q^ 


Qroii  Elec.  Power  (xlO*  MWe)  NBI  Fusion  Power  (xl0~*  MW) 


2.25 
1.75 

1.25 
0.75 

1.7 

1.6 

1.5 

1.4 

1.3 

1.2 


Edge  Safety  Factor,  q^ 


Edge  Safety  Factor,  9.^ 


Ill 


Edge  Sa£st7  Factor, 


Edge  Sa£stjr  Factor,  q^ 


Edge  Sa£et7  Factor,  q^ 


Edge  Saf^  Factor,  9^ 


nccirc.  Power  FVec.  (xlO"’)  TP  Coil  lleiiitive  Power  (MW) 


2 


Edge  SefeQr  Faetaz,  9^ 


Edge  Safety  Factor, 


^3.  125 

S 

I 

I3.O75 


<1 

£3.025 


2.975 


Edge  Safety  Factor, 


W*u  Load  (xlO-'MW/m»)  Diveilor  Heat  Ued  (MW/m*) 


114 


Edge  Sa£sty  Factor, 


1.10  - 1 - 1 - 1 - p 

1.05  - 

1.00  >  - . .  . 

S  0 .95  - 

^0.90  L - 1 - 1 - 1 - L 

3  5  7 

Edge  Sa&t7  Factor, 


Edge  Saiet}r  Factor,  q^ 


Edge  Safstj  Factor,  7^ 


LITERATURE  CITED 


[1]  J.P  Holdren  et  al.  ’’Exploring  the  Competitive  Potential  of  Magnetic  Fusion 
Energy:  The  Interaction  of  Economics  with  Safety  and  Environmental 
Chjiracteristics” .  Puaion  Technology^  13:7,  1988. 

[2]  R.W.  Conn  et  al.  ’’Economic,  Safety  and  Environmental  Prospects  of  Fusion 
Reactors”.  Nuclear  Fusion,  30:1919,  1990. 

[3]  A.E.Dabiri.  ”An  Overview  of  D-®He  Fusion  Reactors”.  Nucl.  Instr.  and 
Meth.,  A271:71,  1988. 

[4]  W.  Kembichler  et  al.  ”D-^He  Fuel  Cycles  for  Neutron  Lean  Reactors”. 

Fusion  Technology,  15:1142,  1989. 

[5]  G.A.  Emmert  et  al.  ” Possibilities  for  Breakeven  and  Ignition  of  D-®He  Fusion 
Fuel  in  a  Near  Term  Tokamak”.  Nuclear  Fusion,  29:1427,  1989. 

[6]  J.F.  Santarius  et  al.  ’’Energy  Conversion  Options  for  ARIES-III  -  A 
Conceptual  D-*He  Tokamak  Reactor”.  In  IEEE  ISth  Symposium  on  Fusion 
Engineering,  Knoxville,  TN,  October  1989. 

[7]  J.R.  McNally,  Jr.  ” Physics  of  Fusion  Fuel  Cycles”.  Nuclear 
Technology /Fiision,  2:9,  1982. 

[8]  L.J.  Wittenberg,  J.F.  Santarius,  and  G.L.  Kulcinski.  ’’Lunar  Source  of  ®He  for 
Commercial  Fusion  Power”.  Fusion  Technology,  10:167,  1986. 

[9]  M.  Kikuchi.  ”  Steady  State  Tokamak  Reactor  Based  on  the  Bootstrap 
Current”.  Nuclear  Fusion,  30:265,  1990. 

[10]  R.S.  Devoto  et  al.  "Projections  for  a  Steady-State  Tokameik  Reactor  Based 
on  the  International  Thermonuclear  Experimental  Reactor”.  Fusion 
Technology,  19:251,  1991. 

[11]  B.B.  Kadomtsev,  F.S.  Troyon,  and  M.L.  Watkins.  "Tolamaks”.  Nuclear 
Fusion,  30:1675,  1990. 

[12]  M.S.  Chance,  S.C.  Jardin,  and  T.H.  Stix.  "Ballooning  Mode  Stability  of 
Bean-Shaped  Cross  Sections  for  High  Beta  Tokamak  Plasmas”.  Physical 
Review  Letters,  51:1963,  1983. 

[13]  Y-K.M.  Peng.  Spherical  Torus,  Compact  Fusion  at  Low  Field, 
ORNL/FEDC-84/7-  Technical  Report,  Oak  Ridge  National  Laboratory,  1985. 


115 


[14]  J.D.  Galambos,  D.J.  Strickler,  Y-K.M.  Peng,  and  R.L.  Reid.  ’’Plasma 
Elongation  Studies  for  ITER- Like  Tokamaks”.  Fusion  Technology,  15:483, 
1989. 


116 


[15]  D.J.  Strickler  et  al.  Equilibrium  Modeling  of  the  TFCX  Poloidal  Field  Coil 
System,  ORNL/FEDC-83/10.  Technical  Report,  Oak  Ridge  National 
Laboratory,  1984. 

[16]  Y-K.M.  Peng  and  D.J.  Strickler.  Features  of  Spherical  Torus  Plasmas, 
ORNL/FEDC-85/6.  Technical  Report,  Oak  Ridge  National  Laboratory,  1985. 

[17]  Y-K.M.  Peng  et  al.  ’’Spherical  Torus:  An  Approach  to  Compact  Fusion  at 
Low  Field  -  Initial  Ignition  Assessments”.  Fusion  Technology,  8:338,  1985. 

[18]  Y-K.M.  Peng  and  J.B.  Hicks.  "Engineering  Feasibility  of  Tight  Aspect  Ratio 
Tokamak  (Spherical  Torus)  Reactors”.  In  16th  Symposium  on  Fusion 
Technology,  London,  September  1990. 

[19]  J.D.  Galambos  and  Y-K.M.  Peng.  "Ignition  and  Bum  Criteria  for  D-®He 
Tokamak  and  Spherical  Torus  Reactors”.  Fusion  Technology,  19:31,  1991. 

[20]  R.T.C.  Smith  et  al.  "START  (Small  Tight  Aspect  Ratio  Tokamak)”.  In  16th 
Symposium  on  Fusion  Technology,  London,  September  1990. 

[21]  A.Sykes  et  al.  "First  Results  from  the  START  Experiment”.  1991.  Draft 
Report. 

[22]  "Quantifying  the  Concepts  and  the  Potential  of  the  Steady-State  Spherical 
Tokamak  (TST)  and  the  Neutron  Spherical  Tokamak  (NST)”.  1990.  Oak 
Ridge  National  Laboratory  Proposal. 

[23]  R.L.  Reid  et  al.  ETR/ITER  Systems  Code,  ORNL/FEDC-87/7.  Technical 
Report,  Oak  Ridge  National  Laboratory,  1988. 

[24]  G.W.  Shuy,  A.E.  Dabiri,  and  H.  Gurol.  "Conceptual  Design  of  a 
Deuterium- ^He  Fueled  Tandem  Mirror  Reactor  Satellite/Breeder  System”. 
Fusion  Technology,  9:459,  1986. 

[25]  S.  Thompson.  "Systems  Code  Cost  Accounting”,  FEDC-M-88-SE-004.  1988. 
Oak  Ridge  National  Laboratory  Memorandum. 

[26]  R.L.  Crane,  K.E.  Hillstrom,  and  M.  Minkoff.  "Solution  of  the  General 
Nonlinear  Programming  Problem  with  Subroutine  VMCON”,  ANL-80-64- 
Technical  Report,  Argonne  National  Laboratory,  1980. 


[27]  Private  communication  with  J.D.  Galambos  and  Y-K.M.  Peng,  Oak  Ridge 
National  Laboratory. 


117 


[28]  0.  Mitarai  et  al.  "Saddle  Point  Ignition  Condition  for  D-®He  Tokamak  Fusion 
Reactor”.  Fusion  Technology,  19:234,  1991. 

[29]  ARIES-III  Parametric/Systems  Calculations  (Version  8.1,  April  1991) 
obtained  through  private  communication  with  C.  Bathke,  Los  Alamos 
National  Laboratory. 

[30]  R.A.  Kreikowski.  "Progress  in  Commercial  Magnetic  Fusion  Energy  Reactor 
Designs”.  Fusion  Technology,  20:121,  1991. 

[31]  W.  Kramer.  "Critical  Safety  Issues  in  the  Design  of  Fusion  Machines”.  In 
16th  Symposium  on  Fusion  Technology,  London,  September  1990. 

[32]  J.F  Sautarius  et  al.  "A  High  Efficiency  D-®He  Tandem  Mirror  Fusion 
Reactor”.  In  IEEE  12th  Symposium  on  Fusion  Engineering,  Monterey,  CA, 
October  1987. 

[33]  V.K.  Rohatgi  and  T.  Vijayan.  "Technical  Issues  in  Fusion  Reactors  -  A 
Review”.  Fusion  Technology,  16:287,  1989. 

[34]  E.E.  Bloom.  "Structural  Materials  for  Fusion  Reactors”.  Nuclear  Fusion, 
30:1879,  1990. 

[35]  N.A.  Uckan  et  al.  International  Thermonuclear  Experimental  Reactor  Design 
Information  Document,  ITER-TN-PH-0-5.  Technical  Report,  IAEA,  Vienna, 
1990. 

[36]  J.  Johner.  "Thermonuclear  Ignition  in  the  Next-Generation  Tokamaks”. 
Fusion  Technology,  19:515,  1991. 

[37]  R.J.  Bickerton,  J.W.  Conner,  and  J.B.  Taylor.  "Diffusion  Driven  Plasma 
Currents  eind  Bootstrap  Tokamak”.  Nature  of  Physical  Science,  229:110,  1971. 

[38]  K.  Miyamoto.  Plasma  Physics  for  Nuclear  Fusion,  chapter  8.  Massachusetts 
Institute  of  Technology,  revised  english  language  edition,  1989. 

[39]  D.P.  Stotler,  R.J.  Goldston,  and  The  CIT  Tezun.  "Ignition  Probabilities  for 
Compact  Ignition  Tolramak  Design  Points”.  Fusion  Technology,  20:7,  1991. 

[40]  B.  Coppi  and  L.E.  Sugiyama.  Questions  in  Advanced  Fuel  Fusion,  Report 
PTP-88/6.  Technical  Report,  Massachusetts  Institute  of  Technology,  1988. 

[41]  G.H.  Miley,  H.  Towner,  and  N.  Ivich.  Nucl.  Eng.  Report  COOO-2218-17. 
Technical  Report,  U.  of  lUinois,  1974. 


